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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. ‘The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in | . Authors are requested to submit MSS. in duplicate. These should be 

using double spacing and leaving a in of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 

in. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zeits. f. Astrophys., 15, 239, 1938. 
A. S. Eddington, Internal Constitution of the Stars, Cambridge, p. 182, Table 24, 1926. 


3. Notation.—Authors should conform closely to the recommendations of Commission 
3 of the International Astronomical Union (Trans. J.A.U., Vol. VI, p. 345, 1938). Council 
has decided to adopt the I.A.U. 4-letter abbreviations for constellations where contraction 
is desirable (Vol. IV, p. 221, 1932). 


4. Diagrams.—These should be drawn about twice the size required in print 
and pr ed for direct photographic reproduction except for the lettering, which 
should be inserted in pencil. Legends should be given in the manuscript 
indicating where in the text the figure should appear. Blocks are retained by the 


Society for 10 years; unless the author requires them before the end of this period they 
are then destroyed. 


5. Tables.—These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts.—When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MS. 


Reading of Papers at Meetings 
8. When submitting papers authors are requested to indicate whether they will be 
willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


9. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1951 APRIL 13 
Professor H. Dingle, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 

Cecil Gilbert, 1 North House, Longbenton, Newcastle-on-Tyne (proposed by 
K. Mitchell) ; 

Bernardo Grossling-Freudenburg, M.Sc., 26 Emperor’s Gate, London, 
S.W.7 (proposed by A. T. Price); 

William Halliday, B.Eng., 3 Hollycroft Avenue, London, N.W.3 (proposed by 
P. F. Everitt); 

Avril Barbara Hart, University Observatory, Oxford (proposed by H. H. 
Plaskett) ; 

John Warren Loveridge, Middleheath, 11 West Heath Road, London, N.W.3 
(proposed by A. Hunter); and 

George Tait, 16 Avenue Road, St. Albans, Hertfordshire (proposed by 
R. W. B. Pearse). 


The election by the Council of the following Junior Members was duly 
confirmed :— 


Bernard Ephraim Julius Pagel, The Observatories, Madingley Road, Cam- 
bridge (proposed by R. H. Garstang); and 
Dennis William Scaima, The Cavendish Laboratory, Free School Lane, 
Cambridge (proposed by R. H. Garstang). 
One hundred and fourteen presents were announced as having been received 
since the last meeting, including :— 
J. B. Sidgwick, Introducing Astronomy (presented by the author); and 
W. M. Smart, The Origin of the Earth (presented by the author). 
The President announced that the Council had elected the following Associates 
of the Society :— 
Dr Gerard P. Kuiper, Yerkes Observatory, Williams Bay, Wisconsin, U.S.A.; 
and 
Professor Polydore Swings, Institut d’Astrophysique, Université de Liége, 
Sclessin, Belgium. 
Professor Antonie Pannekoek delivered the George Darwin Lecture (see 
p. 347), taking as his subject “The Origin of Astronomy”. 
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After the lecture, the President delivered a short address, and presented the 
Gold Medal of the Society to Professor Pannekoek, to whom it had been awarded 
for his work on astrophysics and on the structure of the galactic system. ‘The 
President said :— 


‘*T suppose what strikes us first about the subject of this lecture which we have 
just heard is its remoteness from the kind of paper which we have usually heard at 
our meetings. Yet Professor Pannekoek has dealt with the actual foundations on 
which our more advanced work is built. He has mentioned that in Babylonian 
times they would be prepared to bet five to one on the success of an eclipse pre- 
diction. Anyone who did the same to-day with some of the predictions of modern 
theories would be regarded as a gambler! ‘The debt which the Greeks owed to 
Babylonian observations forms an interesting study, and papers on this subject 
are so few that we are fortunate to have had one at this meeting. 

‘* Professor Pannekoek, it is now my pleasure to hand to you the Gold Medal of 
the Society, which is the highest honour the Society has it in its power to bestow. 
Will you accept it with our thanks for your past work and our best wishes for further 
enjoyment in the many fields of astronomy in which you have gathered so rich a 
harvest.” 


MEE'TING OF 1951 MAY 11 
Professor H. Dingle, President, in the Chair 
‘The election by the Council of the following Fellows was duly contirmed: 
*Richard Myer Baum, Stocks Lane Post Office, 1 Dee Banks, Boughton, 
Chester (proposed by P. A. Moore); 

Ronald Winston Cable, 77 Holland Street, Wembley, Perth, Western 
Australia (proposed by R. W. Boggis) ; 

Selwyn Rhys Davies, B.Sc., 12 Brohanddgar, Felinfoel, Llanelly, South Wales 
(proposed by B. E. Featherstone) ; 

Thomas Charles Mansfield, M.A., 7 Knole Road, Sevenoaks, Kent (proposed 
by H. Spencer Jones) ; 

Tadeusz Boleslaw Slebarski, M.A., University Observatory, St. Andrews, 
Scotland (proposed by E. Finlay-Freundlich) ; 

Stanislaus Alexander Vincenz, B.Sc., A.R.C.S., D.I.C., Imperial College of 
Science, South Kensington, $.W.7 (proposed by J. McG. Bruckshaw) ; 
and 

Andrew Young, B.Sc., Department of Applied Mathematics, ‘The University, 
Liverpool (proposed by E. Finlay-Freundlich). 


One hundred and ten presents were announced as having been received since 
the last meeting, including :— 


J. W. Mitchell, Fundamental Mechanisms of Photographic Sensitivity (presented 
by the editor). 


* Transferred from Junior Membership. 








THE ORIGIN OF ASTRONOMY 


George Darwin Lecture, delivered by Professor Antonie Pannekoek on 


1951 April 13 


‘The origin of astronomy, as of all science, lies far back in prehistoric times. 
Primitive man in his struggle for life could not do without a certain knowledge 
of the surrounding natural phenomena; the better he was acquainted with them, 
the better his life was secured. ‘Thus, when he reached the cultural stage which 
ethnology calls the higher stage of barbarism, a first, merely practical, knowledge 
of physics, chemistry and biology was contained in his techniques, his agriculture 
and his cattle raising. At this stage some knowledge of celestial phenomena 
must also have been present, just as to-day we find it with primitive tribes. Whence 
did it come? Sometimes the opinion has been expressed that it was the beauty 
of the starry heavens, the mysterious silent course of its luminaries with their 
continuous change of aspect, that attracted the attention of man and awakened the 
desire to find out their cause. This opinion is clearly a reflex of modern conditions, 
where physics and chemistry are closely tied up with economic life, whereas 
astronomy stands outside the direct necessities of life and is studied for its own 
sake, for the beauty of its object and the wide scope of its problems. In opposition 
to this Schiaparelli (x) rightly argued that the first study of astronomy was imposed 
upon man by the desire to satisfy the necessities and comforts of his life. Indeed, 
economic necessities were at the root of his astronomical as well as of his physical 
and chemical knowledge. The heavens above him belonged as an essential 
part to the surrounding world which he had to know. 

So, to understand the origin of astronomy, we have to see what practical 
necessities directed the attention of primitive man towards the celestial phenomena. 
‘here is first the need for orientation or direction finding. When commerce 
had developed out of nomadic cattle-herding or coastal fishery, the caravans 
in the desert and the ships at sea needed the celestial luminaries to guide them. 
During daytime the Sun and at night the stars directed their travels. ‘Thus the 
\rabs from olden times were acquainted with a number of stars, especially with 
the 27 groups which the Moon successively visits in its monthly course. In 
Homer’s Odyssey the goddess Kalypso instructs Odysseus in sailing home to 
keep the Bear at his left hand. ‘The Polynesian island dwellers in the Pacific, 
before their contact with Europeans, were well acquainted with the stars; as 
skilled mariners they steered their vessels by means of the rising and setting points 
of different stars, as a kind of celestial compass; they had schools for navigation 
where, by means of globes, young people were instructed in astronomy. 

More important still and more. universal than for this orientation in space 
was the use of the heavenly bodies for what may be called the orientation in t:me, 
the time reckoning, the calendar. All work of man, at least before the industrial 
revolution, was bound up with the alternation of the seasons. ‘The natural 
phenomena on which his activities depend were regulated by the period of the 
solar year. In Northern and moderate climates they depended on the alternation 
of a dead winter and a life-abounding summer, in Southern and tropical countries 
mostly on the alternating of a dry and a wet season. As hunters or fishermen 
they had to follow the migrations of the animals; as herdsmen or farmers they had 


to regulate their work according to the seasonal life-periods of animals and plants. 
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Thus, in the planning of his work, in order not to be misguided by irregularities 
of the weather, man had to have some independent means for counting the days 
and recognizing the right time for sowing and harvesting. Such means were 
provided by the Sun itself, and by other celestial phenomena. 

The phenomena connected with the yearly variation of the Sun’s declination 
are eo ipso exact indications of the seasons, and were recognized and used as such. 
Many primitive peoples to-day are found to use the solar phenomena as signs of 
the time of year. Dayak tribes measure the shadow of a vertical stick held at 
arm’s length to know the month (2). The Eskimos and the Zuni Indians are 
reported to observe the solar solstices by watching the moments of extreme 
rising and setting at points of the horizon which are marked by stones (3); such 
moments are important dates celebrated by religious rites. 

Another important characteristic of the season is given by the stars, most 
manifestly by their heliacal rising, in the morning just before sunrise. The 
Australian aborigines at Torres Strait for the beginning of sowing await the 
appearance of a bright star called Kek, probably Canopus or Achernar, for which 
a close watch is kept (4). Analogous practices are mentioned for other primitive 
peoples. ‘That the same is true of prehistoric man is indicated by fragments of 
the oldest historical! testimonies. Everybody knows how in ancient Egypt the 
heliacal rising of Sirius in the month of June was taken to announce the flooding 
of the Nile. From ancient Greece, Hesiod’s poem Works and Days mentions 
for the different rural activities what phenomena of the stars are connected with 
them; e. g. “ When the Atlas-born Pleiades rise, then begin mowing, but ploughing 
when they disappear ”’ (5). 

But more than the Sun and the stars, the Moon occupied the chief place in 
time-reckoning. Its regular changes of phase offered the most practical period 
for counting the days. Everywhere, with all peoples—even those who after- 
wards adopted, as we do, pure solar reckoning—the Moon’s synodic period of 
294 days was originally the basis of the calendar. The month begins with the 
first appearance of the new sickle Moon as a slender arc in the Western evening 
sky and as a consequence the day begins with the evening. ‘The Moon’s gradual 
waxing and waning impressed man as the life-history of a living being, ever 
again new-born. In the sacred books of many peoples it is said that the Moon 
has been created for the purpose of measuring the time. The Moon was the 
oldest and most venerable of the gods; its first appearance was attentively watched 
and hailed, and still more its fullness, when it dispels the darkness of the night, 
was celebrated with religious ceremonies. 

Thus a certain knowledge of celestial phenomena, as a part of surrounding 
nature, is present in prehistoric man, involved with all his other technical know- 
ledge in the practical activities of sustaining his life. Often it was a speciality 
in the hands of old people or priests, who perfermed the small functions of 
spiritual leadership in these communities. We cannot truly call it astronomical 
science, any more than his techniques of making tools and preparing food is called 
physical or chemical science. Science in the real sense of the word could only 
arise when man entered into a higher stage of development, the stage of 
civilization. 

The transition from the prehistoric stage of barbarism to civilization is marked 
by the invention of script. This means more than that from now on written 
history supplants the oral tradition of legends. ‘The importance of script is 
that it gives visible form to the conceptions that before had only a spiritual 
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existence in our mind. ‘The words of the spoken language expressing these 
conceptions had already been in use for countless generations, as a means of 
mutual intercourse and incentive to action, but only as part of instinctive 
practical life; once the fleeting sound had passed, nothing material remained. 
Through script, however, conceptions acquire an existence of their own, visible 
and permanent; they can be handled and stored, compared and connected with 
one another. Now science as a system of generalized knowledge, of relations 
between general notions and abstract conceptions, becomes possible. Science 
is theory; it deals with things in general, with the abstractions formed by our 
mind out of the phenomena. Thus the rise of science as theoretical knowledge 
must be seen as part of the larger historical process, the rise of civilization. 

The origin of civilization has to be sought in the realm of social developments. 
‘‘War and Class’’, writes the English historian Arnold Toynbee, ‘are the two 
congenital diseases of civilization” (6). Now war was inherent in the preceding 
stages too, during savagery and barbarism. But ‘‘Class’’, i.e. the separation 
of a ruling class from the working masses, has been the social characteristic of 
civilization only, in all its forms. Its role in the development of astronomical 
science must now be traced. The countries to which we have therefore to turn 
our attention are Babylonia and Greece. 

The first civilizations arose about 5000 years ago in the fertile plains of Egypt, 
Mesopotamia, India and China. Here the silt deposits of the great rivers, Nile, 
Euphrates, Hoang-ho, offered a good soil of marvellous fertility for a numerous 
population, but continual care was needed to regulate the water supply by dykes 
and canals. At regular or irregular moments the river broke through the dykes, 
flooded the land and deposited new layers of silt. This care could not be left 
to the separate villages or districts with their often opposed interests and so a 
strong central power was necessary to provide for the common interest. The 
multitude of small communities with their own chieftains and deities consolidated 
into large monarchies. ‘This was rendered possible because the fertility of the 
soil procured abundant surplus produce for the sustenance of a separate class 
of ruling officials. 

Moreover a strong state power was needed for defence against the poor 
warlike tribes which in the surrounding mountains and deserts could find a 
meagre living only. They made it their business to attack and plunder their 
prosperous neighbours. So a division of labour was necessary; an army of 
soldiers had to protect the farmers and in time became a ruling group, with their 
chief as king taking into his hands the judicial and administrative powers. 
Sometimes the plundering tribes became conquerors and settled as a ruling 
aristocracy among the subjected farmers; theirs was now the task of staving off 
new aggressors. ‘The result was the same in both cases. Beside the military 
power stood the civil officials. With the separate townships organized into a 
large empire their local deities too were organized into a pantheon and their local 
priests into a powerful hierarchy, the spiritual leaders of society. 

With the appearance of these ruling classes new needs arose. Their greater 
wealth expressed itself in luxury, in the development of refined arts, in spiritual 
culture. The conditions were now present for the introduction and spreading 
of script. The ruling monarchs wanted to glorify their deeds, their victories 
in war and their temple-building in peace-time. Missives with reports and 
instructions were interchanged between the king and his officials and governors, 
and his decrees were fixed as laws by writing them down. Writing as a special 
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capacity of priests and scribes was also put at the service of private individuals 
for their business contracts. Thousands of such contracts as well as other 
documents, written on clay in cuneiform characters, have been dug up from 
temple ruins of the Babylonian towns and are preserved in our museums. 

What was the effect of the new conditions upon the rise of astronomy ? 
A special group of officials had now to pay attention to the celestial pheno- 
mena. ‘The regulation of time was one of the chief tasks of the priesthood; from 
their staged towers they had to observe the first appearance of the new sickle Moon 
and to announce the new month. Then the Moon reckoning had to be adapted to 
the seasons of the year. Because 12 periods of the Moon are 11 days shorter than 
the solar year, every 3 or 2 years a 13th month must be intercalated, not for 
any theoretical or scientific reason, but in order that the offering of the first crops 
to the gods may take place at the prescribed time, at the full Moon ceremonies of 
the first month Nisannu. In such primitive societies the important agricultural 
activities were at the same time religious festivals; economics and worship, 
politics and religion formed an indissoluble unity. The service of the gods was 
bound to strict ritual; the calendar, the time regulation of worship, was a sacred 
duty of the priests, the source of their prestige and spiritual power. Originally 
the intercalation of a 13th month took place according to practical necessity ; 
when at the end of the 12th month it appeared that ripeness of the crops and the 
beginning of harvesting could not be expected in the next fortnight, the 12th month 
had to be repeated. ‘The inscriptions show that about 2000 B.c. the intercalations 
took place at irregular intervals, and hence were probably determined by practice. 

This may have caused annoyance, but a better way presented itself. When 
looking out for the new Moon the priests must have perceived that in every 
succeeding month the surrounding stars had changed, gradually moving farther 
West till they disappeared, being replaced by more easterly constellations. 
At the same time new stars rose in the morning twilight. So the phenomena of 
the stars indicated the seasons and could afford a more regular and exact means 
for intercalation of the additional months. There are many texts where the names 
of constellations are connected with the names of months, indicating their use 
for calendar purposes. An easy rule is given in a text published by George 
Smith: ‘‘ When at the first of Nisannu the Moon and star Mulmul (the Pleiades) 
stand together the year is common; when at the third day of Nisannu the Moon 
and the star Mulmul stand together the year will be full” (7). The latter case 
means that the first month falls so early in the seasonal year that addition of a 
13th month becomes necessary. As a consequence we find that long before 
1000 B.c. the Babylonians were acquainted with a considerable number of 
constellations, by names in inscriptions and by pictures on boundary stones. 
They are to a great extent identical with our own, evidently transmitted by the 
Greeks. ‘That in these same centuries also the brightest of the planets, with its 
alternation as evening and morning star, was perceived, is evident from the famous 
so-called Nindar-anna tablet from about 1600B.c., where its phenomena are 
described, and their significance. 

For in the course of time it was no longer only, or chiefly, calendar needs that 
determined man’s interest in the stars. Ever more astrology, the doctrine of 
the significance of the celestial phenomena for the events on Earth, occupies his 
mind. Now what happens in the sky is followed with keener interest as signs 
of the destiny allotted to people and states. Belief in signs and presages is as 
universal in primitive man as his belief in invisible spiritual powers surrounding 





ATTN io 


No. 4, 1951 The origin of astronomy 351 


him and influencing his work and life. ‘To win their favour, to find out their 
intentions, to avert or appease their hostility, by magic charms and incantations 
is an important part of his daily practice. Most of these spirits have their abode 
in the heavens ; and especially in these Mesopotamian plains, with their resplendent 
sky, the stars impressed man’s mind by their radiant brightness. Thus in early 
times the conception of a close connection between the stars and man’s fortunes 
had already risen in the mind of the Babylonian priests. Old Sumerian inscrip- 
tions speak of the “favourable star”’ for temple building. This belief was alive 
through all the centuries when Babylon, as a rich commercial town and the 
capital of a large empire, had become a centre of culture for the surrounding 
world. But it intensified to an all-pervading fervour when in the centuries 
after 1000B.c. Assyria rose and became the most powerful state in the Near 
East, extending its dominance not only over Babylon and all Mesopotamia but 
also over Syria, Palestine and even Egypt. In the vicissitudes of world politics 
at the time and in the variable fortune of war the need was strongly felt for aids in 
foreseeing the future. So the astrologers in direct service of the king had to give 
their omens before every great enterprise, and from all the great sanctuaries the 
priests had to send regularly their reports on what they had seen in the heavens 
and the meaning thereof. From the ruins of the library of King Assurbanipal 
ten thousands of more or less damaged sherds with cuneiform inscriptions have 
been dug up, many of them containing astrological reports. ‘These enable us to 
obtain a clear picture of the character and scope of their astronomical observations. 

The chief objects observed were the Moon and the planets, because by their 
variable aspect and irregular wanderings they offered the greatest variety of 
phenomena suitable for astrological interpretation. ‘The Moon was now observed 
with more attention in all its phases; especially were its risings and settings 
observed in the middle of the month, about full Moon. When the Moon is full 
on the night of the 14th, the normal time, it was a lucky omen; when full Moon 
happened on the night of the 13th, 15th or 16th, it was abnormal, hence a bad 
omen. Here astrology and calendar were merged; deviation in the calendar 
was considered an unlucky sign and had to be restored at the end of the month. 
The eclipses, of course, were important phenomena, mostly unfavourable; but 
as to the different sides of the Moon that were successively obscured and the 
different hours of the night, they had significance for different countries, for 
Akkad, the southern Babylonia proper, for Elam in the East, for Amurru in the 
West. So the details were often so well noted that nine centuries afterwards 
Ptolemy could use them in the construction of his lunar theory. 

But the most important effect of astrology was that now besides the Moon 
the planets strongly captured the attention of man. They were seen wandering 
between the stars in the most intricate and unexpected ways, as living beings 
roaming through the starry landscape. ‘They are the stars of the great gods 
who rule the world and manifest themselves in these bright luminaries. Their 
appearances and disappearances, their motions through the constellations of the 
zodiac, their stationary points and retrogradings, their meeting with one another and 
with bright stars, presented an almost endless variety of phenomena. Although 
from earlier times we have already a text with omens from Venus, it is in these 
later centuries of Assyrian power that a high tide of astrology provides us with 
an abundance of reports containing omens from all the five planets. There we 
read (8): ‘“‘ Mercury went back as far as the Pleiades’’ ; “‘ Jupiter enters Cancer ”’ ; 
‘** Venus appears in the East”; ‘‘ Mars is very bright’’; “Jupiter appears in the 
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region of Orion’’ ; ‘‘ Mars stands in Scorpio, turns and goes forth with diminished 
brilliancy”’; ‘‘ Saturn has appeared in the Lion’’; “‘ Mars approached Jupiter” ; 
and soon. There is not a trace of scientific interest in these texts; the mind of 
the reporters is entirely occupied by the omens: when such or such happens, 
“it is lucky for the king, my lord”’; or, ‘‘ copious floods will come” ; ‘there will 
be devastation”’; ‘the crops will be diminished”; “the king will be besieged ”’ ; 
“the enemy will be slain’; “there will be raging of lions and wolves”; “the 
gods intend Akkad for happiness’’; and soon. Yet, with all those observations, 
these reports represent a considerable astronomical activity. For the first time 
in history a large number of data on the planets had been collected; it implies a 
detailed knowledge of facts about their motion. 

This does not mean that we have here the beginning of a science of astronomy. 
Mere facts do not form a science. Facts are the basis of science. Science is 
the systematization of facts into general conceptions and rules. Now it might 
be imagined that these ancient observers, by comparing their observations with 
old records and looking for regularities, must gradually have detected them. 
But we must beware of regarding these Babylonian priests as astronomers like 
ourselves, animated by the same spirit of scientific purpose and scientific research. 
Such a state of mind certainly was entirely absent then. It is the very problem 
of the origin of science to discover how it originally came about that, without 
knowing about scientific aims, people were able to build up scientific theory. 

The circumstance that made this possible for astronomy was the occurrence 
of extremely simple and striking periodicities in the celestial phenomena. What 
looked irregular on occasional and superficial observing revealed its regularity 
in a continuous abundance of data. Regularities were not sought for; but 
regularities imposed themselves, without giving surprise. They aroused certain 
expectations. Expectation is the first unconscious form of generalized knowledge, 
like all technical knowledge in daily life growing out of practical experience. 
Then gradually the expectation develops into prediction, an indication that the 
rule, the regularity, has entered consciousness. In the celestial phenomena the 
regularities appear as fixed periods, after which the same aspects return. Know- 
ledge of the periods was the first form of astronomical theory. 

Did the astrologers in Assyrian times already possess such knowledge? 
There are some texts indicating that they did. In one of them (g) we read: 
“To the king, my lord, I sent: an eclipse will take place. Now it has not passed, 
it has taken place. In its happening the eclipse portends peace for the king, 
my lord.” Here it is stated that a former prediction has been confirmed. The 
prediction of a lunar eclipse indeed is very simple, because 5 or 6 such eclipses 
follow one another always regularly with 6 months’ interval; and when this 
series has ended a new series starts 11 or 17 months later. Hence, if an eclipse 
has been observed we may bet five to one that 6 months later another will take 
place. There is a difficulty, however; on the average half of them will be 
invisible because they occur in the daytime, when the full Moon is below the 
horizon. This irregularity must have hampered the easy discovery of the 
continuous series by the Babylonian priests. But their own regular observations 
of the Moon’s risings and settings in the midst of the month must have shown them 
a way towards understanding such gaps. When full Moon happens in the 
night, then in the morning the Moon will set only after sunrise; Sun and Moon, 
god Shamash and god Sin, are seen at the same time, one East, the other West. 
If however, the Moon does not wait for the Sun but sets before sunrise, full 
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Moon has still to come; and if in such a case an eclipse was expected, it will 
remain invisible. That the Babylonian priests were aware of these conditions is 
shown by a published text, probably from the 7th century B.c.: ‘The eclipse 
passes, it does not take place. If the king should ask: what omens hast thou seen ? 
—the gods have not been seen one with another” (10). Here clearly the eclipse 
was predicted, and now the reason is given why it must remain invisible. 

In the same way some regularities may have been perceived with the planets, 
and may have resulted in a certain expectation. Thus in an Assyrian text we 
read: ‘Jupiter has stood for a month over its reckoned time” (11). Hence its 
time had been reckoned. The most regular phenomena, such as the periodical 
invisibility of the planets in the vicinity of the Sun, or, for Jupiter and Saturn, 
their stately course along the zodiac, making a circumference in 12 and in 29 years 
respectively, must have become practical experience for the astrologers, just as 
the sequence of the Moon’s figures is practical experience for the savage. We may 
be sure of it without having express statements in written records. Some 
centuries later we find texts expressly testifying to knowledge of planetary periods, 
already at a higher stage. But, in the meantime, conditions had changed 
considerably. 

First politically, when those warring little countries of old had been absorbed 
into the inner peace of the huge Persian empire, the old omens of lucky for Akkad, 
unlucky for Elam, lost their sense. ‘Then astronomically, the detection and 
gradually arising knowledge of the periodicity and the periods must have brought 
about a new attitude of these observers towards the celestial phenomena. ‘They 
no longer were the ignorant earthlings anxiously looking at the sky for the messages 
which the supreme masters of the world should inscribe there. ‘They knew in 
advance what would be inscribed there, they could predict it, they had some 
kind of spiritual power over it. Their capacity of prediction gave them a new 
social qualification, as people who know the ways of the gods. Thus their 
observing activity impelled them to greater zeal and persistence as a ritual duty 
in the service of the gods who reveal themselves in these stars. It became more 
conscious and complete, more detailed and precise; distances of the planets to 
bright stars were noted numerically, hence perhaps measured, although we know 
nothing of the instruments used. 

The result of this now genuinely astronomical practice appears in the rather 
scanty inscriptions of the following centuries of rule by the Persian and the 
Seleucid kings. They show that now the larger periods of the planets were 
known, the common multiples of synodic period and period of revolution: 8 years 
for Venus, 71 and 83 for Jupiter, 46 years for Mercury, 47 and 79 years for Mars, 
59 years for Saturn; after these intervals the phenomena repeat themselves at 
nearly the same places. These periods moreover were practically applied for 
prediction purposes, in such a way that the phenomena, observed the appropriate 
number of years earlier, were transferred, with small corrections, to the year in 
question. ‘Thus there exists a tablet entitled: ‘‘'The first day, the phenomena, 
the motions and the eclipses for the year 140”’ (i.e. Seleucid Era) (12); it gives 
phenomena of Jupiter for the years 69 and 57, of Venus for 132, of Mercury for 94, 
of Saturn for 81, of Mars for 61 and 93: exactly the years in which the data will 
produce the phenomena for the desired year by addition of the periods named. 
Ephemerides computed in this way have been found for different years, up to 
the beginning of our era. Here we have astronomical theory in the form of 
knowledge of the periods consciously applied in prediction. 
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A still higher and more perfect form of theoretical knowledge was reached by 
Babylonian astronomy in the last centuries B.c. Although Babylon had lost its 
greatness and the Chaldaean country, dominated by Parthian kings, was cut off 
from the new centres of commerce in the Mediterranean, astronomical activity 
persisted and now reached its highest stage of perfection before it was extin- 
guished. ‘The tablets dug up from the ruins of old city temples of this time, 
mostly much-damaged fragments, show nothing but rows of numbers arranged in 
columns, intermingled with the names of months and zodiacal signs. Their 
structure and meaning have been unravelled by the work of Father F. X. Kugler, 
S.J., both an able astronomer and an Assyriologist. In the planetary tables each 
column contains the successive values of the longitude and the date of one of the 
special phenomena, opposition positions and heliacal rising and setting. All 
qualitative description, such as proximity to a star or entrance into a constellation, 
is absent; there is only pure quantity, given with utmost formal precision in 
sexagesimal fractions. Now for the first time, because all the data are given by 
mere numerical values, the inequalities in the planetary motions can become 
evident as a periodical increase and decrease of the successive intervals. ‘The 
Chaldaean astronomers did not represent them, as we do, by sinusoidal wave 
functions based on geometrical construction, but by purely arithmetical zig-zag 
functions, running with constant differences up and down between an upper and 
a lower boundary where they are reflected. ‘The same procedure was applied 
in their lunar tables, where time and place of the conjunctions and oppositions, 
the first sickle Moon and the eclipses were computed by means of exact data on 
the Moon’s different periods. 

We are confronted here not only with far more detailed knowledge of the 
planetary and lunar motions and their periodic inequalities than in the previous 
records, but also with its presentation in the form of more abstract scientific 
theory. ‘The previous form was an excellent device for technical prediction 
based on the knowledge of periods; but for every prediction reference to former 
observations was necessary. Here the tables stand by themselves, as pure theory, 
no longer needing any observational datum. ‘They have not the meaning of mere 
predictions; the rows of numbers embrace past and future undistinguished and 
can be extended in both directions indefinitely. It is abstract theory in the form 
of concrete data; instead of by rules stated in words it is given by the regularities 
in the tables of figures. ‘The curious fact is that we do not know how it originated ; 
these tables appear suddenly, without intermediate forms that could tell us about 
their relation to former phases of knowledge. This, however, is only part of 
the general problem of how the different stages of Babylonian astronomy have 
developed out of or beside one another and from what observational data, a problem 
hardly tackled as yet. 

Babylonian astronomy in the last centuries presents an admirable system of 
organized, i.e. theoretical, knowledge. But it is merely formal theory, entirely 
devoid of any physical interpretation. It is numerical knowledge only; nowhere 
is the question of structure raised. Its mathematics is only arithmetic, with no 
geometry. It appears that for these Chaldaean astronomers the planets were not 
solid bodies moving in orbits through space, but rather, as in olden times, luminaries 
wandering—although now with unravelled regularity—along the firmament. 
‘The cause of this restricted character must be sought for in the fact that this 
science was entirely in the hands of priests, who by a mighty power of tradition 
were bound to the old forms of worship and ritual. All their astronomical work 
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was service to the gods: the most modern planetary tables we find to be introduced 
by the sacral formula: “ In the names of god Bel and goddess Beltis, my masters, 
an omen.’’ So the way to further theoretical progress, to structural physical 
theory, was barred. If we consider structural theory an essential characteristic 
of theoretical science, then for the origin of astronomy one step more must be 
made. For this we have to turn our attention to another people, to the Greeks. 

In living conditions and in character Greece and the Greeks, country and 
people, in every respect present the greatest possible contrast to Babylonia. 
In this mountainous peninsula, deeply indented by the sea, with good harbours 
but with only small bits of arable land separated by wild-wooded ridges, the 
people remained divided in small communities with their local deities and local 
priests. ‘They became seafarers and merchants, and many of their expanding 
numbers emigrated and settled on foreign shores, in Asia Minor and Southern 
[taly. In such traders and colonists a freer, more independent and more daring 
spirit developed than among farmers staying at home under permanent conditions ; 
they were less bound to tradition, more open to new ideas. Commerce gave 
rise to industry, the fabrication of wares for export, and inventiveness became 
the basis of industrial progress. Here, in the Greek colonies, there arose a ruling 
class of free and wealthy citizens, the first men in history comparable to the 
burgher class that in modern times placed its stamp upon European society. 
Here a new spiritual life revealed itself in a flowering of poetry and of new philo- 
sophical world-conceptions. 

These first Greek philosophers amaze us as much by the boldness of their 
ideas about world structure as by the scantiness of their knowledge of astronomical 
phenomena. From Xenophanes and from Heraclitus—-famous even now as 
che first thinker who saw the world as a process of eternal flux——the opinion is 
reported that the Sun and the stars were formed anew every day. ‘The cosmical 
theories of these times were not devised to explain celestial phenomena but 
their object, as shown by B. Farrington (13), was to apply known technical 
processes to the world at large. ‘These Greek philosophers were thinkers, not 
observers. Simple phenomena such as the identity of evening and morning 
star or the course of the planets through the zodiac, which with a little attention 
could have been noticed by anyone, are attributed as great discoveries to famous 
philosophers. And this situation remains throughout the centuries in which 
Greece reached the summit of political power and of literary and artistic brilliancy. 
Plato, who knows and describes the planets in their sequence, their colour and 
brightness, says of their wanderings that they are incalculable in multitude and 
marvellously intricate. Eudoxus* famous theory of the homocentric spheres 
was an attempt to explain the mere qualitative fact of retrograde motions, which 
he did by means of a most ingeniously devised system of uniform rotations. 
But it did not fit the retrogressions as they factually were ; obviously there were no 
observational data to test the details of the theory. 

Judging from the amount of their observations and detailed knowledge of the 
celestial luminaries, the Greeks in their most glorious times cannot be said to 
have been great in astronomy. But they did something else. ‘They developed 
geometry and a geometrical mode of thinking as a most admirable system of 
demonstrable truths, fit to become the basis for a future higher stage of astronomy. 
‘They developed a conception of the world, or rather a new way of looking at it, 
as a system of material bodies and spheres moving and revolving in three- 
dimensional space. They stood within this world as artificers dominating it 
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mentally, looking from the outside at the constituting bodies with the Earth 
itself in the midst, handling them as objects of geometry, cutting their spheres 
with planes and making the circles and triangles objects of geometrical demon- 
stration. In Aristarchus’ determination of the relative distance of Sun and 
Moon we have a fine example of this attitude towards the heavenly world, Greek 
astronomy at its best. Yet in the work expounding this method the apparent 
diameter of Sun and Moon is assumed to be 2°, four times too large. 

A new epoch set in when as a consequence of Alexander’s conquests the 
Greek and the Oriental world were merged, the Hellenistic era. Now the Greeks 
got access to the observing practice and to the abundance of observational knowledge 
of the Babylonians. ‘The Oriental example gives a new impulse; observations 
are made at Alexandria and elsewhere. But what we know of it deals only with 
the fixed stars, not with the planets. And Ptolemy afterwards states that 
Hipparchus from lack of a sufficient number of observations could not bring his 
theory of the planets to completion. So far, Greek astronomy is seen to keep 
its former character, but now the Greek thinkers had at their disposal the numerical 
results of Babylonian astronomy. Now their geometrical world schemes could 
be filled in by exact data on periods and inequalities ; now the numerical perfection 
of Chaldaean science could assume the form of geometrical structure in space. 
The result was the epicycle theory, the supreme product of ancient astronomy. 

In the epicycle theory we may be said to have reached our goal, true astro- 
nomical science. There is a conception of science which considers only the 
latest truth as real, true science. It sees the epicycle theory as a primitive 
erroneous system, to be superseded 1700 years later by the true world system. 
Looking, however, from an historical point of view we see that it represents 
correctly the relative motions of the planets in their circular orbits, and that only 
the original assumptions as to the zero point of the motions had to be corrected 
afterwards. It is scientific theory in the strictest sense of the word, systematiz- 
ation of observational facts in a world structure suited to the computation and 
prediction of future events. The study of the history of science presents it to our 
view as a most important stage in the growth of knowledge, conditioned by the 
entire social development of the time, as the outcome of observing and thinking 
activities during a long series of centuries, and as the starting point of later progress 
towards modern science. Thus we understand how at so early a time, when in 
no other realm of natural phenomena had knowledge risen above the stage of 
technical experience, astronomy could rise to the high rank of a science. 
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RADIO EMISSION FROM THE ANDROMEDA NEBULA 
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Summary 


A detailed account is given of the measurement of the radio-frequency 
radiation from the Great Nebula in Andromeda (M31) using a wave-length 
of 1°89 m. It is shown that the measured intensity is in good agreement 
with that to be expected on the assumption that the generation of the radiation 
is similar to that in the Galaxy. From these results it is estimated that the 
total extra-galactic radio emission amounts to about 1 per cent of the whole 
radio flux observed on the Earth. 





1. Introduction.—The experiments of Jansky (1) in 1931 showed that radio- 
frequency radiation is reaching the Earth from the general direction of the Milky 
Way. Subsequent measurements by Hey, Parsons and Phillips (2), by Reber (3) 
and by Bolton and Westfold (4) have shown that the intensity contours of this 
radiation correspond well with the structure of the Galaxy derived from visual 
observations. Although the source of the radiation remains unknown, it is 
generally accepted that the major part of the radiation must be generated by some 


mechanism which is widespread in the Galaxy. 

Attempts have been made in the past to discover whether or not the extra- 
galactic nebulae generate similar radio emissions. For example, Reber (3) used 
a paraboloid of 30 feet aperture on a wave-length of 1-87 metres, but failed to find 
any conclusive evidence of radiation from the nebula M31. During the last 
few years, however, a much larger paraboloid of aperture 218 feet has been 
available at the Jodrell Bank Experimental Station. Calculations of the 
intensity to be expected from M31 based on the assumption that it radiates in 
a similar way to the Galaxy suggested that with this paraboloid, used in conjunction 
with the best available receiver, it might be possible to detect the radiation if 
it existed. 

This paper describes the results of an experiment designed to test if such 
radio-frequency radiation is being emitted by the Great Nebula in Andromeda 
(M31). A short account of the results has been published previously (5). 

2. Description of apparatus.—Fig. 1 shows a block diagram of the apparatus. 
The aerial is a paraboloid of diameter 218 feet and focal length 126 feet. The 
reflecting surface is formed by long wires which are spaced 8 inches apart and which 
run parallel to one fixed plane of polarization. 

Electromagnetic waves of wave-length large compared with 8 inches and whose 
polarization is parallel to the wires are reflected by this surface and are focused 
on to the primary feed which consists of two dipoles and their reflectors mounted 
on a central mast 126 feet high. This mast is pivoted at its base in such a way 
that it may be tilted in the north-south plane to 15 degrees on either side of the 
vertical. The latitude of the system is N.53° 14’ and it is therefore possible to 
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direct the beam to declinations between N. 38° and N.68°. ‘The direction of the 
beam is estimated by measuring the angle through which the primary feed is 
displaced from the vertical and applying to this angle a correction established from 
experiments (6) carried out with small paraboloids. 


» 
La 
< « ws 
! NOISE 
MATCHING 
~ UNIT 1 GENERATOR 
pa yd c 
' 


SELECT VE DWASE 
ed SENSITIVE 


















































RECORDER 5 * RECEIVER 
FREQUENCY 











AMPLIFIER RECTIFIER 


7 


















































a 
sS 


Block diagram of apparatus. 


The gain and beam-width of the aerial have been calculated for various 
wave-lengths, taking into account the polar diagram of the primary feed, the actual 
shape of the reflecting surface and the reflection coefficient of the system of 
wires. ‘They have also been measured at 4-2 metres using an aircraft as a trans- 
mitting source. ‘The beam-width at 1-89 metres has been measured by obsery- 
ations of the radio source in Cygnus. ‘The values obtained are shown in 
Table I. 

"TABLE | 
Characteristics of Aerial System 
72 Mc./s. 158-5 Me.’s. 
‘Theoretical Measured Theoretical Measure: 
Gain (over half-wave dipole) 690 700 2130 aoa 
Beam-width to half power 4° 30° 4° 20° 2’ 00’ 1° 58 


Fig. 2 shows the shape of the beam measured in both planes by observations 
of the source in Cygnus, that is with the beam tilted about 14 degrees from the 
vertical. ‘The broken line in Fig. 2(a) shows the theoretical beam shape when 
the beam is directed to the zenith. 

The aerial system is connected to the receiver through 294 feet of coaxial 
feeder which has a loss of 2-7 decibels at 1-89 metres. 

The receiving equipment is based on that described by Ryle and Vonberg (7). 
Their arrangement has been modified to allow the measurement of aerial 
temperatures below room temperature. A rotating switch S, (Fig. 1) revolves 
at about 1200 r.p.m. and switches the receiver alternately to a noise generator A 
and to the aerial in parallel with a noise generator B. Both noise generators 
are tungsten filament diodes type CV 172. Any difference between the power 
received in the two positions of the switch appears as a low-frequency modulation 
of the receiver output. ‘This low-frequency component is amplified by a selective 
low-frequency amplifier and is applied to a phase-sensitive rectifier operated from 
a switch S, which is synchronous with the aerial switch. The output of the phase 
sensitive rectifier controls the power delivered from the noise generator B which 
is in parallel with the aerial. ‘The system automatically maintains a balance 
such that the power from the aerial plus that from the generator B balances the 
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power from the generator A. ‘The power from the generator A is kept constant 
and the variations in the power from the aerial are recorded as variations in 
power output from the generator B. ‘These variations are recorded continuously 
by a standard recording milliammeter on a moving chart. 

In the region of sky around M31 the effective aerial temperature remained 
below room temperature and, except when calibrating, the output from the 
generator A was simply that of a resistance at room temperature. ‘I'he load 
impedance of the generator A is accurately pre-set and an adjustable matching 
unit is used to equalize the impedances of the two branches at the switch 5). 
As a further precaution the output of the generator B is calibrated by means of 
the generator A. 














(a) (b) 
. 2.-~Polar diagrams of 218-ft. aperture paraboloid. (A= 1°89 metres.) 
(a) Beam shape in right ascension. 
(6) Beam shape in declination. 
Ordinates: Intensity in arbitrary units. 
Abscissae : Degrees. 

The receiver has a midband frequency of 158-5 Mc./s. determined by a 
crystal-controlled local oscillator. ‘The pre-detector band-width is 1-5 Mc./s., 
and the radio-frequency amplifier is a “cascode” (8) with a noise factor of 3°5. 
The time constant of response of the whole system is about 10 seconds which 
is small compared with the 10 minutes taken by a point source at declination 
N. 40° to transit the aerial beam. 

The minimum power incident on the aerial, which can be detected by such 
an equipment, can be estimated theoretically, assuming that the limit is set by 
the random fluctuations inherent in the equipment, and not by external inter- 
ference. 

Assuming that 

N =noise factor of the receiver, 
B =band-width of receiver in cycles per second, 
t =time constant of response of the equipment, 
T =total effective temperature at input to receiver, 
AT =r.m.s. fluctuations of output meter, 


at~nr, |(Z). 


it can be shown (g) that 
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Substituting the parameters of the equipment described above, AT ~ 0-7 deg. C. 
Practical experience shows that under the conditions of this experiment the 
minimum deflection which can be detected on a single trace is about twice the 
r.m.s. fluctuation of the output indicator. Using the theoretical value of AT 
derived above, the minimum change of aerial temperature which can be detected 
should be about 1-4 deg. C. 

Allowing for the loss in the feeder cable and for the decrease in gain of the 
aerial system when the beam is tilted 14 degrees from the zenith, this value of AT 
corresponds (at 1-89 metres wave-length) to an incident randomly polarized flux 
of 10-*5 watts/square metre/c.p.s. Practical experience shows that the minimum 
detectable flux depends upon the time of day and on the weather. During the 
day, and in fine weather, the noise level of the equipment is controlled by reception 
of the “‘ quiet’’ Sun in minor lobes of the aerial, and by man-made static. Charged 
rain, electrical disturbances in the atmosphere and radiation from the “disturbed” 
Sun overload the equipment. During the night the noise level falls off gradually 
until about midnight and increases again at sunrise. Between midnight and 
sunrise the minimum detectable flux appears to be controlled sometimes by the 
theoretical limits discussed in equation (1) and sometimes by unidentified signals 
which appear to vary with the weather and which probably have a terrestrial 
origin. On favourable nights it is about 10~*> watts/square metre/c.p.s., which 
is the limit estimated theoretically. By averaging the results of several inde- 
pendent observations it is possible to extend this limit by an amount which has 
not yet been established by a quantitative experiment. In the present work the 
limit has been successfully reduced to 5 x 10-*° watts/square metre/c.p.s. by 
averaging four records. 

3. Method of observation.—The experiment was carried out by fixing the 
beam at a number of different elevations corresponding to the region around 
M31, so that for each diurnal rotation of the Earth the beam swept out a strip of 
sky 2 degrees wide in declination and 24 hours in right ascension. For each 
elevation of the beam the intensity of the radio flux received was recorded continu- 
ously for about 100 to 200 hours. Altogether, recordings were made for 
several elevations of the beam, corresponding to declinations between N. 38° 47’ 
and N. 43°00’. The variation with right ascension of the flux received at each 
declination was found by averaging the results obtained on several successive 
sweeps. 

4. Results—The total number of records obtained by the method described 
above was 90, of which 20 were spoiled by the effects of heavy rain which persisted 
throughout 1950 August and September. Of the remaining records 50 showed 
the existence of a localized source of radio-frequency radiation at 00" 40 R.A. 
Plate 4 shows one of these records taken with the beam at declination N. 40° 11’. 
The remaining 20 were either taken when the beam was directed to declinations 
N. 38° and N. 43°, or interpretation of the records was made difficult by inter- 
ference in the period 23"30™R.A. to or®40™R.A. The curves shown in Fig. 3 
were constructed from 30 records which were undisturbed by interference. 
In order to show clearly the change in intensity with right ascension along each 
sweep the records have been displaced by an arbitrary amount, corresponding to 
different declinations, so that the ordinate represents only the relative intensities 
for the points in each sweep. The sweeps between declinations N.39° 48’ 
and N. 41°58’ indicate a localized source of radiation about 00" 40™ R.A., super- 
imposed on a slow change of intensity which reaches a minimum between 
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Facsimile of a single record of the transit of the source in Andromeda taken on 1950 October 31 
with the aerial beam directed to declination 40°11'N. 
Ordinates : Intensity in arbitrary units. 
Abscissae : Hours and minutes of right ascension. 
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00" 20™ and 00830™R.A. This slow change corresponds to a gradient of the 
galactic background radiation, which reaches a minimum when the beam is near 
its maximum southern galactic latitude. 


43° 00’ 
Being 41° 58’ 


Gethin a en 
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Fic. 3.—Jntensity recorded during transit of the source in Andromeda with the aerial beam 

directed to different declinations. 











Ordinates : Relative intensity in arbitrary units. 
Abscissae : Hours and minutes of right ascension. 


These curves cannot be converted into true contours of the intensity from this 
region because the change of aerial impedance with beam tilt did not permit the 
difference in absolute intensity between observations at different declinations to 
be measured reliably. ‘Therefore the background intensity was extrapolated to 
oo" 40™ R.A. for each curve and the intensity of the background radiation at this 
point was taken as the zero for each curve when plotting the contours shown in 
Fig. 4. The contour system (Fig. 4) therefore represents the intensity relative 
to the background intensity at 00" 40™R.A., the gradient in declination of the 
galactic radiation having been removed. The gradient in right ascension has 
not been removed and so causes a distortion of the contours. 

Analysis of the results yields the data shown in Table Il. The source has a 
finite apparent width in right ascension and declination, as can be seen by comparing 
the apparent widths with the corresponding widths of the intense source in 
Cygnus, which is at approximately the same declination as the observed source, 
and whose apparent width is known to be less than i’-5 of arc (10). This 
comparison eliminates any uncertainty due to possible distortion of the beam. 
The appropriate curves are shown at (a) and (c) in Fig. 5. Record (5) in Fig. 5 
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shows the intensity received from another weak source at approximately the same 
declination as the observed source. 
The intensity given in Table II is an integrated value over the source and is 


twice the intensity observed on one polarization, as it is assumed that the radiation 
a4 


39 


38 
1-10 C100 00-50 00 40 00 30 00 20 0010 
Fic. 4.—Contours of radio-frequency flux observed near the source in Andromeda with a 2-degree 
beam. (1 unit=10~** watts/square metre/c.p.s. A=1°89 metres.) 

The contours do not represent the absolute intensity of the radio flux at each point. As explained 
in the text the gradient of background flux in declination has been removed and the contours show 
the intensity for each point above the background flux at oo4 40 R.A. The gradient'of background 
flux in right ascension has not been removed and therefore distorts the contours. 

The broken line shows the outline of the nebula derived from a photograph. 

Ordinates : Declination (degrees north). 

Abscissae : Right ascension (1 division represents 10 minutes). 


Taste II 
Data on Source in Andromeda 

Right Ascension 005 40™ 158+ 308 
Declination 40° 50’ N.+ 20’ 
Apparent widths to half-power 3° (dec.) x 34° (R.A.) 
Apparent widths to half-power of the source 

in Cygnus 2}° (dec.) x 2° (R.A.) 
Intensity (integrated over source) 10~*4 watts/sq. metre/c.p.s. + 25 per cent 

















2020 940 1920 
Fic. 5.—JIntensity recorded during transit of three sources through the aerial beam. The scale 

of intensity for each diagram has been adjusted so that the widths in right ascension can be 
compared easily. 

(a) Source in Andromeda. 

(6) Weak source in Perseus. 

(c) Intense source in Cygnus. 

Ordinates : Intensity (arbitrary units). 

Abscissae : Hours and minutes of right ascension, 
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is randomly polarized. The previously published (§) value of 4 x 10~*5watts/square 
metre/c.p.s. was calculated assuming the apparent diameter of the source to be 
small compared with the beam-width, and radiation from the outer regions of the 
source was neglected. Further measurements and analysis showed the extension 
of the source to be so great that this assumption was not justified. 

The intensity of the source is so small compared with the background 
radiation and with the random fluctuations in the receiver that the difficulty in 
interpreting the records is the main cause of errors in the datain Table II. There 
is no significant systematic error in right ascension as the position of the beam 
in this coordinate was checked by reference to the intense source in Cygnus 
whose right ascension is well established. As the declination of the source in 
Cygnus is not yet known with satisfactory accuracy (10) it was not possible to 
check the position of the beam in declination in the same way. This introduces 
an additional error into the determination of declination, as the relationship 
between the angle of tilt of the mast and the angle of the beam is not yet known to 
better than I per cent. Therefore there may be a systematic error of the order 
of +10’ in the measurements of declination. The error in measuring the actual 
tilt of the mast was negligible. 

The error in the measurement of intensity is caused partly by the weakness 
of the source and partly by inaccuracy in the calibration, these two errors being 
of comparable magnitude. 

5. Discussion of results.—The coordinates of the observed source lie close to 
the centre of the Great Nebula in Andromeda M31 (R.A. 00" 40™, declination 
N. 41°00’, Epoch 1950). The apparent dimensions of the source are consistent 
with a source of radio-frequency radiation of size comparable with the main body 
of the nebula M31.* It is possible that these results could be due to a fortuitous 
grouping of two or more point sources with their effective centre coincident with 
that of M31. From a consideration of the number of observed point sources of 
intensity comparable to 10~*4 watts/square metre/c.p.s. and, assuming a random 
distribution of sources, it can be shown that the probability of observing two 
sources simultaneously in the beam is about ;4,. It is therefore improbable 
that the finite width of the source is due to a coincidence of two point sources. 

The celestial coordinates and finite size of the source indicate with a high 
degree of probability that the source may be identified with the nebula M3r. 
It is assumed that the two companion nebulae NGC 205 and NGC 221 (which 
have apparent photographic magnitudes of 10-8 and 9:5 respectively compared 
with the apparent photographic magnitude 5 of M31 (11)) are not contributing 
substantially to the radiation received. 


* Note added in proof.—The apparent extension of the source can, in principle, be deduced by 
applying to the observed results a correction for the finite beam-width of the serial. However, 
the low intensity of the source in comparison with the gradient of background radiation precludes 
an accurate analysis of the extension, particularly in right ascension. ‘The values for the extension 
of the source given in a previous publication (§) were underestimated. A subsequent analysis 
showed that the observed variation of intensity with declination could be accounted for by a source 
in which the intensity is a maximum at the centre and decreases to 10 per cent of the central intensity 
at points separated by 150’. The apparent extension of the source in right ascension appears to be 
of the same order, but the gradient of background radiation makes this result unreliable. 

Inspection of the contour system given in Fig. 4 shows that, when allowance is made for the 
distortion in right ascension of the contours by the gradient of background radiation, it is possible 
to interpret the results as due to an elliptical source orientated in a similar manner to M31. 

The authors are obliged to Mr J. H. Piddington of the Radiophysics Laboratory, Sydney, whose 
criticisms have led to the inclusion of this footnote. 


25* 
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6. Comparison of the radio emission from M31 and the Galaxy.—From the 
results described above it appears evident that radio emissions from M31 are 
being received on the Earth. In this section we compare the measured intensity 
of this emission with that to be expected if M31 is similar to the Galaxy. An 
estimate of the intensity to be expected from M31 has been made in two ways, 
(a) by consideration of the luminosities and (b) by analysis of the contours of 
radio flux from the Galaxy. 

(a) Comparison of radio flux and light flux.—This method makes the assump- 
tion that the ratio of the radio flux from M31 to the radio flux from a volume near 
the Sun equals the ratio of the light flux from M31 to the light flux from the 
same volume near the Sun. 

The aggregate luminosity per cubic parsec near the Sun is known to be 
0-045 suns. ‘The radio flux from the same volume can be estimated from the 
temperatures observed (4) in the directions of the galactic poles, if it is assumed 
that the radio flux originates in bodies distributed in a manner similar to the 
visible stars and that absorption is negligible. ‘The temperatures of the galactic 
poles have been used because in these directions the estimates of stellar distri- 
bution are the most reliable. 

Consider an aerial with a beam of solid angle w, effective area A and gain G 
directed towards a galactic pole and receiving radiation of wave-length A metres. 

Let 

p,=the density of sources per cubic metre at a distance r metres from the 
Sun, 
Po=the density of sources near the Sun, 

47a =the total power radiated by an average source (watts/c.p.s.), 

I =the intensity of the radio flux at the Earth (watts/square metre/c.p.s.), 
P =the power received by the aerial (watts/c.p.s.), 
T =the effective temperature of the aerial (deg. K.), 


a 
4m 
where R is the distance from the Sun to the boundary of the Galaxy in a direction 


perpendicular to the galactic plane, the boundary being taken as the distance at 
which the stellar density falls to 1 per cent of that near the Sun. Therefore 


rR 
P=kT= af | i p, adr, 


“R 
then [= |, are 


rR 
4nr? dr = aw | Pedr, (2) 


where k is Boltzmann’s constant (the factor 4 is included because the aerial 
receives radiation in one polarization only), but 


_47A _4n 
G= =<: (3) 
Hence 


9 
= 


kT= 


R 
2 - p, ar. (4) 


The integral | p, ar, evaluated by using the star density/distance tables given 
0 


by Bok (12), equals 1-4 x 107 po. 
Taking (4) T =600 deg. K. and A=3 metres, 
&py = I°3 x 10-“ watts/steradian/cu. metre/c.p.s. 
= 3°8 x 10° watts/steradian/cu. parsec/c.p.s., 
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it has been shown that the law connecting the intensity (J) of the background 
radiation with wave-length (A) is J =constant x A*, where 0-35 <a<0°65 (10). 
By assuming J: :A°®, 
%&py = 3 X 10° watts/steradian/cu. parsec/c.p.s. 
at a wave-length of 1-89 metres. This is the radio flux per steradian from a cubic 
parsec near the Sun. 

The ratio of the light flux from M31 to that from a cubic parsec near the Sun 
is 2:2 x10! (The luminosity of M31 is taken as 10° suns and that of a cubic 
parsec near the Sun as 0-045 suns.) 

Assuming this same ratio for radio-frequency radiation, then the total radio- 
frequency radiation from M31 is 2-2ap, x 10’ watts/steradian/c.p.s. The 
intensity of this radiation at the Earth should therefore be 

2'2%p, x 101° 
[2 
where / is the distance of M31 in metres. 

If 1/=230,000 parsecs, then J\y,, = 1-2 x 10~*4 watts/square metre/c.p.s. which 
is of the same order as the observed value. 

(b) Analysis of contours of radio flux from the Galaxy.—The contours of radio- 
frequency radiation over most of the celestial sphere have been measured by 
Bolton and Westfold (4) at a wave-length of 3 metres. These have been used to 
construct a rough model of the Galaxy as it would appear from the distance of M31. 

If absorption is negligible, then the intensity of the radiation observed from 
this distance will be approximately independent of the orientation of the Galaxy 
with respect to the observer. In this analysis the Galaxy was assumed to be 
viewed along the equatorial plane. The observed contours were projected on 
to a cross-section of the Galaxy through the galactic centre and perpendicular 
to the observer. The Galaxy was assumed to be a spheroid with dimensions 
30,000 x 5,000 parsecs and to have a distribution of mass symmetrical about its 
minor axis. The projection was made on the assumption that the radiation from 
any volume of the Galaxy is proportional to the mass in that volume, From these 
projected contours {7(Q)dQ was evaluated over the whole cross-section of the 
Galaxy, where 7(Q) is the effective temperature observed over an elementary 
solid angle dQ subtended at an observer at a distance of 230,000 parsecs. The 
intensity (J,,) observed at this distance will be * 


watts/square metre/c.p.s., 


2kr 
Ig= | T(Q) dQ. 


* This expression can be derived as follows :— 


Consider a beam of solid angle w where w is large compared with the apparent diameter of the 
source. The effective temperature (7,) of this beam will be given by 


T= f7(Q) dQ 
w 


using the same notation as in the text. The power (P) received by the beam is given by 
RIT(Q)d2Q Ig A 
w apa 


P=kT,;= 


i 2kf7(Q) dQ 
56 _— 
From equation (3) Aw=A’, 
sp 
*. Ig= x | 72) dQ. 
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Evaluating the integral as described above, the value of J, is found to be 
g x 10-“ watts/square metre/c.p.s. at a wave-length of 3m. Then, assuming 
I::°5 as before, the intensity on the wave-length of 1-89 m. used in the present 
work should be 7 x 10-* watts/square metre/c.p.s. for the Galaxy observed at 
the distance of M31. 

By assuming, as above, that the radio flux is proportional to the mass, then the 
intensity to be expected from M31 can be estimated by comparing its mass with 
that of the Galaxy. Estimates of the mass of M31 vary from 3 x 10 to 10" solar 
masses compared with an estimate of 10!" to 2 x 10" for the Galaxy. The intensity 
of the radio-frequency radiation from M31 observed on the Earth would thus be 
expected to be within the range I x 10-** to 7 x 10-*4 watts/square metre/c.p.s. 
which, in fact, includes the observed value. 

The agreement between the theoretical values calculated above and the 
measured intensity suggests that as far as radio-frequency emission is concerned 
M3I possesses similar characteristics to the Galaxy. 

7. The total radiation from extra-galactic nebulae.—The contribution to the 
radio flux incident on the Earth from extra-galactic nebulae has been calculated 
for a very simple model of the universe, assuming that the ratio of radio flux to 
light flux found for M3r is true for all nebulae independent of type. 

The intensity of the radio flux from M31 observed on the Earth at a wave- 
length of 1-89 metres is 10-* watts/square metre/c.p.s. and therefore assuming it to 
be an isotropic radiator the total radio flux emitted is 5 x 10!® watts/steradian/c.p.s. 
M31 has an absolute visual magnitude of — 17-5 compared with — 15-2 for the 


average nebula. Therefore the radio flux from an average nebula will be about 
5 x 1018 watts/steradian/c.p.s. 


If | o@=the average space density of nebulae (assumed isotropic), 
4my =the total radio flux from an average nebula, 
I, =the intensity at the Earth of extra-galactic radiation per steradian, 
then if radiation is received from nebulae in a sphere of radius D and the effects 
of absorption and recession are neglected, 


Ip=yoD (from equation (2)). 


Putting y =5 x 101° watts/steradian/c.p.s., o = I-g x 10-!* nebulae per cubic light- 
year, then for a universe of radius 5 x 10° light-years (corresponding to the range 
of the 100-inch telescope), [= 4°75 x 10-*4 watts/square metre/steradian/c.p.s. 
If the radius of the visible universe is taken as 1-8 x 10° light-years, as suggested 
by the apparent recession of the nebulae, 


I, =1-7 x 10~-** watts/square metre/steradian/c.p.s. 


This latter value of J, corresponds to an effective aerial temperature at 
1-89 metres of about 2deg.K. ‘The average aerial temperature observed over 
the whole celestial sphere at this wave-length is about 200deg.K., and this 
calculation therefore suggests that the extra-galactic radiation contributes about 
I per cent of the total radio flux incident on the Earth. 
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Summary 


Summaries are given of the atomic ionization and recombination data 
required in the calculation of interstellar ionization equilibrium. The 
chemical composition of the interstellar gas is discussed from the standpoint 
of seeking to establish whether there is any evidence for discrepancies 
between the interstellar and ‘mean cosmic” abundance ratios. The 
important Ca/Na ratio is the only case for which there is an apparent 
discrepancy, but it is shown that this may be removed if recent work on the 
ultra-violet radiation field is accepted as correct. 





1. Introduction.—In order to obtain information on the chemical composition 
of the interstellar gas it is necessary to study the ionization equilibrium which 
is established. Strémgren (1) has shown that the procedure appropriate to 
the conditions in interstellar space is to equate the rate of photo-ionizations from 
the ground state of a given ion to the sum of the recombination rates to all 
excited states. Both of the rate coefficients involved depend upon the continuous 
atomic absorption cross-sections, and in view of recent work on these it appears 
to be desirable to tabulate the atomic data in a form which is most convenient 
for astrophysical applications, and which is independent of any particular 
assumptions about the conditions prevailing in interstellar space. 

We denote by Ry: and Ry. the number density of ions of a given element 
X in two consecutive stages of ionization, and by R, the number density of 
electrons. The rate of recombinations per unit volume is then given by 


gh =aRyR, (x) 
For an electron temperature 7, the radiative recombination coefficient « is 
given by 
I 2\2 w*, . r oO . 
a= (=) (mkT ME a exp (J AT.) | 0 ay(hvPexp(—Av/kT,)d(hv), (2) 
where wi, I* and a* denote the statistical weight, ionization potential and 
continuous absorption cross-section of the Ath excited state of X’ and wy» is 
the ground-state statistical weight of X”. 
The rate of photo-ionizations per unit volume is given in terms of the energy 
density of radiation p, and the ground-state absorption cross-section a} by the 
relation 


c 
g!™ = Ry. fat 5 Pedy = Ry T. (3) 


We refer to [' as the ionization coefficient. Equating the ionization and 


recombination rates we obtain 


R,.Ry TT 
a ort (4) 


which is the Strémgren ionization equation. 
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2. Atomic data for the calculation of ionization and recombination rates 

2.1. The recombination coefficient.—The recombination coefficients for a 
number of ions of importance in interstellar space are given in Table I. In 
calculating the sum of partial recombination coefficients, all states for which no 
entry has been made in the table have been taken to be hydrogenic. This 
approximation has been discussed by Bates, Buckingham, Massey and 
Unwin (2); the accuracy of the present atomic data will be discussed in 
Section 2.3 below. For interpolating the temperature variation of « the most 
convenient procedure is to plot «7} against log T,,. 

2.2. The ionization coefficient.—Estimates of the interstellar radiation field 
are usually expressed as a sum of black-body energy densities multiplied by 
corresponding dilution factors, 


p= & 8(Ti)p,(T;) (5) 


It is therefore convenient to tabulate the atomic data in the form of black-body 
ionization rates given by 


y(T;) = 1-00 x 108 exp (—1/kT;) F(T;) sec, (6) 


with F(T,) =10 | a\(I +)? exp(—e/kT;) de, (7) 
0 


where a! is in cm.”, J and « in Rydbergs (13°5, e.V.), 7; in deg. K. and & in 
Rydbergs/deg. K. (i.e. k =0°638 x 10-°). The interstellar ionization rate is then 
given by 

P= 2 8(7;)(T;) sec. (8) 


The function F(7;) is given for a number of ions in Table II. In certain 
cases it is necessary to consider the ionization rate with a cut-off for energies 
greater than the ionization potential of hydrogen ((J+e)>1); calculations 
with a cut-off have been marked with an asterisk. 

2.3. Accuracy of the atomic data.—The accuracy of most of the atomic data 
should be adequate for treating the problem of interstellar ionization 
equilibrium (9). For the hydrogenic data the correction factor g in Kramers’ 
formula has not been taken fully into account, but the figures given should be 
correct to within about 10 per cent. For the neutral atoms the only case for 
which the use of the hydrogenic recombination coefficients is likely to introduce 
any appreciable error is for Ca I, for which the s, p and d series are markedly 
non-hydrogenic. It is possible that the total « in this case is somewhat 
underestimated. For the positive ions the departure from hydrogenic values 
is more marked; in particular it will be seen that the Ca II 3d recombination 
term is some six times greater than hydrogenic. Some corrections have been 
made for Na II and K II. 

In a number of cases it has been necessary to extrapolate the experimental 
or theoretical cross-section curves. ‘The only case for which this introduces 
some uncertainty is for the ground state of K I. The experimental results of 
Ditchburn, Tunstead and Yates (11) give the cross-section out to I700A. 
((1+«)=0°54 Ryd.), where it is still rising steeply. The uncertainty 
introduced is greatest for the higher temperature ionization rates, more than 
70 per cent of the K I ionization coming from the extrapolated part of the 
curve for 7, greater than 15,000 deg. K. It should be possible to obtain a reliable 
theoretical cross-section at higher energies using the method of Bates (16). 
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TABLE I 
Final Final Recombination coefficient « 
ion state in units 10~"* cm.°/sec. Reference 
T,=10° T,.=10° T,=10° 
HI 0°16 0*50 1°6 (3) 
0020 0077 0°24 (3) 
0°057 0°22 069 (3) 
0°0059 0°027 0°084 8 (3) 
0018 = 0'079—Ss«O8'25 (3) 
0'020 ©0'092—Ss«0'29 (3) 


0°0039 0021 0°065 (4) 
O'OIO) =—«_« 0053 o°17 (4) 
o’oI2 0062 0°20 (4) 
0°0062 0'033.—s«OvlI (4) 


0°095 0°83 4'°9 (I, 2, 5, 6) 
0°41 2°0 8-6 

0°16 0°47 1'°5 (7) 

0°33 1 | 7°6 

00002 00013 00041 (8, 9) 
0°030 =: 0095 0°30 (10) 
©°00004 0°0003 0°0009 (10) 
0009S s«O"0441 o'13 (10) 
0°17 1°2 5°9 

©'00I 0000, (11) 
0'025 0°039 (12) 
0°16 I'l 5°6 

0°031 (13) 
o°17 

0°65 . (3) 
0°087 . (3) 

0°26 ‘ (3) 
0°033 (3) 

o'10 ; (3) 

or12 : (3) 
0°021 (4) 
0056 (4) 
0°064 (4) 
0°034 (4) 

0°74 (1, 2, 5, 6) 
2'2 

0°19 (14) 
0°50 “ See text 
1°6 

0°39 ’ ; (5) 

0°69 . “ See text 
2°0 

0°69 (9, 15) 
0°003 (9, 13, 15) 
0°17 . (15) 

1°7 
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TasBLe II 
Atomic data for the calculation of ionization rates 
F(T) 

TI in Ryd. F(50,000) F(25,000) F(15,000) F(10,000) 
1°000 160 9°0 5°6 3°8 
0828 12°0 9°4 7°0 51 
0°378 0'0081 0°0067 0°0056 0°0046 
0°319 0°42 o'll 0°026 0°007I 
0°449 10°4 6°1 4°0 2°7 
3°47 330°0 150°0 86-0 560 
2°34 640°0 300°0 170°0 120°0 
0°873 0°57 o’21 o'lr 0067 
0873 0°13 0°095 0073 0°055 

For references to atomic data see Table I, 
The asterisk denotes cut-off for (J+) >1. 


3. Interstellar ionization equilibrium.—To make further progress we must 
make a definite choice for the interstellar radiation field. In Table III we 
give the ionization rate [ and the equilibrium constant I'/« calculated using 
the temperatures and dilution factors given by Dunham(r7). In using these 
data there are four possible sources of error to be borne in mind: 

(i) Dunham states that his results are of a preliminary nature and that 
further work may call for some revision. 

(ii) The correction for general galactic absorption and scattering is rather 
uncertain; we have used Dunham’s data without an absorption correction. 

(iii) The ultra-violet radiation field is obtained using assumed black-body 


curves for stars of representative spectral classes. This source of error will be 
discussed in Section 4.3 below. 

(iv) In discussing the conditions in a single cloud, allowance must be made 
for the possibility of marked local variations from a statistical field due to the 
influence of nearby stars. 


Taswe III 
Ionization equilibrium data using the radiation field of Dunhamt 
r D/a 
Ion in 10718 sec,~* in cm.~* 


T.=10° T,=10% T,=10* 
H!I 43 105 22 5°0 
cr 110 320 63 14 
Nal 4°0 24 3°3 0°68 
KIr® 42 260 38 7°5 
Cal 1400 8200 1200 230 
Na Il 0°62 0°39 0°079 0°018 
Kil 27 13 3°0 o'7I 
Ca Il 2°5 1°5 o°31 o’071 
Ca II* 0°86 O's! oll 0°025 

As before, an asterisk denotes a cut-off for (J+ ¢) >1. 


t+ Some of the data in this table may be compared with a similar table given by Strémgren (1). 
‘The only marked discrepancy is for Ca II, for which Strémgren gives values of ['/a about four times 
greater than ours. This is probably due to the use of the cross-section of Bates and Massey (13) 
without correction for a normalization error pointed out by Green (15). 

The opportunity is taken of correcting two errors in a previous paper (§) by the present author, 
Firstly, in Table III, T,=10, 10%, 10‘ should read 7,=10‘, 10°, 10%. Secondly, an incorrect 
multiplying constant has been used in calculating the K II recombination coefficient. Corrected 
values are given in the present paper. 
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4. The chemical composition of the interstellar gas 


4.1. Introduction.—There are not at present sufficient data available to- 
make possible an unambiguous determination of the chemical composition 
of the interstellar gas. It would seem reasonable to expect, however, that the 
composition would be the same, to a first approximation, as that of the Sun, 
stars, meteors and gaseous nebulae, as expressed in tables of ‘‘ mean cosmic 
abundances”. We will therefore adopt the approach of seeking to establish 
whether there is any evidence for systematic departures from the mean cosmic 
abundance ratios. In order to have a definite basis for discussion we adopt 
as standard the cosmic abundances given in a _ recent review by 
Harrison Brown* (18); values for the elements of interest are quoted im 
Table IV. 

TABLE IV 
Mean cosmic abundance ratios (after Harrison Brown) 
Relative number of atoms 
H 10° K 
Cc 230 Ca 
Na E*3 Fe 
Mg 25 Ti 
Si 29 
C+ Mg+ Si+ Fe=336 

Assuming the composition of the interstellar gas to be the same as the meam 
cosmic ratios suggest, the electron density in regions of ionized hydrogen 
(HII regions) may be equated to the hydrogen density, and in regions of 
neutral hydrogen (HI regions) it may be equated to the combined densities 
of C, Mg, Si and Fe. A guide to the temperature is provided by the recent 
work of Spitzer (20), which suggests that the electron temperature is of order 
10* deg. K. in H II regions and 10* deg. K. in HI regions. For the ionization 
rates in H I regions we take, where necessary, a cut-off for (J+«)>1. 

4.2. The sodium density.—We consider first the observational results for the 
density of neutral sodium atoms. For an HI region most of the sodium will 
be singly ionized, and the total sodium density is therefore given by 


I r 
Rya=Ryau=Ryai- ° (F) Ss 


x 


Assuming cosmic abundance ratios, the hydrogen concentration is then given by 


R R r 
8. (7) (2 - 
m (a) (et) rer Gi, I 
=I°OxI0°Ry,; (T,=100 deg. K.). 


From a very extensive analysis of the sodium D absorption lines Spitzer (21) 
deduces the concentration of neutral sodium atoms per cm.* to be 5 x 10-8 in 
the more opaque clouds and 2 x 10-* in the less opaque clouds.f These figures 
give a hydrogen concentration in the clouds of 9/cm.* and 6/cm.*, which is 
in good agreement with the value derived by Strémgren (1) from the extent 

* The problem of the chemical composition of the gaseous nebulae has been further discussed 
by Aller (x9). 

t It is to be noted that the neutral sodium would be expected to be strongly concentrated in 


clouds (proportional to the square of the hydrogen concentration). We therefore consider Spitzer’s. 
data for the clouds rather than that for an average over all space. 
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of hydrogen emission regions. Our discussion of this point follows closely the 
treatment of Strémgren. 

It is to be noted that the derived hydrogen density is proportional to the 
square root of the ionization function, and will therefore be insensitive to the 
approximation used for the radiation field. The resulting electron concentration 
is of order 0:003/cm.*; carbon, which gives most of the free electrons, is readily 
verified to be almost completely ionized under these conditions. The hypothesis 
that the absorbing clouds are HII regions leads to a hydrogen concentration 
of order one atom per cm.* in the clouds, which appears to be too low. 

4.3. The Ca/Na ratio.—Many writers have commented on an apparent 
discrepancy between the interstellar value for the Ca/Na ratio and the mean 
cosmic and terrestrial values.* It is of particular interest to see whether this 
discrepancy is real. It appears to be well established that the interstellar 
Ca II1/Na I ratio is of order unity. Strémgren (1) obtained the value 0-6 as 
a mean for seven stars. In Table V we give results for a further 14 stars 
suitable for analysis by the doublet ratio method (1), and common to lists of 
equivalent widths for NaD lines given by Spitzer (21), and of equivalent 
widths for Ca II H and K given by Spitzer, Epstein and Li Hen (24). 


TABLE V 


Star R cats ; Star R cat = 
HD number R ya1 Weight HD number’ R yay Weight 


1337 1°3 °'9 184915 o’7I 2°0 
35619 0°36 2°0 185780 2°2 1°o 
163472 0°29 O's 186994 0°39 1°o 
164353 o°31 are) 188209 0°29 03 
167264 0°83 1'O 189550 0°40 o'7 
167838 0°56 o's 190919 2°32 orl 
175876 o'7I | are) 204172 1°3 3°0 


Weighting factors have been assigned according to the weights given in the 
equivalent width tables, with an additional factor depending on the part of the 
curve of growth involved. The weighted mean gives Ro, 1;/Ryq1=0°85, in good 
agreement with Strémgren’s value. We adopt the value o-8 for future 
calculations. There appears to be little doubt as to its essential correctness. 

For an H I region with T= 100 deg. K. the ionization data in Table IV give 


Rea o8 1% + wey 


Rya F R, + 0°68 

For R, of the order of 0-003/cm.* this gives Ry,/Ry,=0°033, which is much 
less than the mean cosmic value of 1-5. The result is quite insensitive to the 
exact value of the electron concentration and to the assumed electron 
temperature.t 

There is a much greater sensitivity to the radiation field, since the NaI and 
Ca II ionization rates depend on quite different spectral regions (2413 — 1800 A. 
for Nal and 1044-911. for Call). This may be expressed in terms of the 
temperatures of the stars responsible for the ionizations. For Na I 80 per cent 
of the ionization comes from stars with 7;< 16,000 deg. K., while for 
Ca II only 20 per cent of the ionization comes from these stars. 

* Eddington (22), Struve (23), Dunham (17), Strémgren (1). 


+ For an HII region we obtain Roqg/Rya=0-06 for Ry=1. For Ry=10/cm.' the ratio is 0-25, 
but this gives the Na/H ratio to be too great by a factor of 100. 
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The possible inaccuracies in the atomic data are certainly too small to account 
for the discrepancy (9), and it is therefore necessary either to accept an 
interstellar Ca/Na ratio much less than the mean cosmic value, or to make a 
rather drastic alteration in the assumed radiation field. There is some evidence 
that the second alternative might be justified. Firstly, for the Sun it is possible 
to calculate the Na I ionization rate using ultra-violet intensities obtained from 
rocket spectrograms, and it is found that the result is about one-tenth of that 
obtained from the conventional 6000 deg. K. black-body model (25). Secondly, 
for stars of spectral type around Bo, which give the major contribution to the 
Ca II ionization, it is possible to calculate the ultra-violet flux using the method 
of model atmospheres. The results of Underhill (26) and of Pecker (27) give 
an ionization rate which is some seven times the value obtained using the 
effective temperatures adopted by Dunham. If these changes can be regarded 
as representative, the previously calculated interstellar Ca/Na ratio must be 
multiplied by a factor of order 70, giving Roq/Ry_=2°3, which removes any 
discrepancy between the mean cosmic and interstellar values. 

4.4. The interstellar lines in  x* Orionis.—The spectrum of x” Orionis is 
particularly rich in interstellar absorption lines. In the first column of 
Table VI we give the ionic abundances as derived by Strémgren (1) from data 
of Dunham (17). These are calculated on the assumption that the absorbing 
path-length, D, is 15 parsecs, which is of the order of the diameter of a single 
cloud.* 
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TaBLe VI 
Abundances from x* Orionis absorption lines 





HII region, 


H I region, T,,=100 deg. K. T.me deg, K. 


Ion densities 
D=15 psc. 





D=15 psc. D=25 psc. D=15 psc. 





NaI 2x107* 
Ki 3xse" 
Cal 2x10~° 
Till 5 x10-* 


x” Orionis 


60 x 10~* 
10 X 10~* 
40 X 107% 

5 Xx 107% 


Cosmic 
70 
go X 10~* 
14 X 107° 
130 X 107% 
5 xX107* 


x? Orionis 
60 x 10-* 
10 X 10 


60 x 10-8 
3x107° 


Cosmic 
40 

50x 107° 
8x 10-* 
80 x 107 
sxze0* 


x? Orionis 


10 X 10~* 
4X10 
4xX10-* 
>5 x 10-* 


Cosmic 
6 
8x 107% 
1X 107% 
10 X 107% 
04 X 107% 



































For an HI region we may assume that nearly all the titanium is singly 
ionized (the ionization energies are 0-502 Ryd. for 7iI and 1-005 Ryd. for 
Till). Equating the total 77 density to the observed 7711 density, and assuming 
cosmic abundances, gives Ry =70/cm.? and R,=0-023/cm.*. That these values 
are greater than would be expected for an average cloud is not surprising 
since the x? Orionis lines are abnormally strong. Using an electron density of 
0:023/cm.° the concentrations of Na, K and Ca may be calculated. The results 
given in Table VI are seen to be in reasonable agreement with the mean cosmic 
values (fitted to agree with 7i). For an assumed path-length of 25 parsecs the 
agreement is seen to be almost perfect; the same agreement would be obtained 
using a path-length of 15 parsecs and increasing all ionization rates by a factor 


* The high-dispersion measurements of Adams (28) suggest that the absorption is almost 
entirely due to a single cloud. 
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of 25/15. If the radiation field were changed in the manner suggested in 
Section 4.3, it would still be possible to obtain the 77, K and Na abundances 
in agreement with the cosmic values, provided a much greater path-length was 
assumed, but the Ca/Na ratio would be a good deal greater than the cosmic 
value. It is not necessary, of course, to assume that the field for a single cloud 
is the same as the statistical field for many clouds. 

We have also considered in Table VI the hypothesis that the absorbing 
cloud is an HII region. For this case it is only possible to give a lower limit 
for the 7i density, since there would be an appreciable amount of 7iIII. We 
have taken R,= Ry =6/cm.*, which gives the best fit for Na, K and Ca. It will 
be seen that the agreement with the cosmic values is not so satisfactory as on the 
HI hypothesis.* 

4.5. The Call/Cal ratio.—Observations of the CalIlI/CalI ratio would be 
valuable in leading to estimates of the electron concentration which would not 
involve any assumptions about the absorbing path-length or the chemical 
composition of the gas. Unfortunately, when the CalI lines are observable 
the Call lines are nearly saturated and an accurate determination of the 
number of atoms is very difficult. Struve (23) has estimated the ratio to be 
between 100 and 150 for 55 Cygni and x Orionist, using assumed equivalent widths 
and curves of growth. What appears to be a more reliable estimate for 
x’ Orionis has been made by Dunham (17), who gives the ratio to be 3500. 
For an HI region (T,=100 deg.K.) this gives an electron concentration of 
0:066/cm.°, which in view of the uncertainties involved, does not compare 
unfavourably with our previous estimate of 0-023/cm.3 (D=15 psc.). For an 
HII region the corresponding value is 2-3/cm.*. 

5. Discussion and conclusions.—It appears to be well established that the 
relative abundances of the elements in the interstellar gas are of at least the same 
approximate order of magnitude as for other cosmic sources. The agreement 
is rather better on the assumption that the average cloud is an HI region than 
that it is an HII region. The only case for which there is some evidence of a 
departure from cosmic values is for the Ca/Na ratio, but in this case it appears 
at least plausible that a revised estimate of the ultra-violet radiation field might 
clear up the discrepancy. Our results apply most directly to the metals, Na, 
K, Ca and 71, and indirectly to the elements, mostly C, which supply the free 
electrons. Hydrogen densities derived from other considerations are in good 
agreement with values derived on the assumption that the composition is in 
agreement with mean cosmic values. 

In discussing the chemical composition of interstellar matter it is necessary 
to take into account the composition of the grains as well as that of the gas. 
According to van de Hulst (29), the grains constitute 0-6 per cent of the toral 
mass of interstellar matter, while for cosmic abundances the metals contribute 
0-3 per cent of the total mass. Appreciable departures from cosmic abundances 
in the gas might therefore be expected if there were a marked concentration 
of metals in the grains. However, recent work{ suggests that dielectric grains, 

* Strémgren (x) has considered the possibility of the Na, K and Ca lines originating in an 
HII region, and the 7% lines originating in an HI region. However, it would appear that an 


HI region of sufficient extent to account for the TiII lines would also be sufficient to account for 
all the other lines. 


¢ Bates and Massey (13) have considered the theoretical implications of this. 
tH. C. van de Hulst (29), Greenstein (30), Whitford (31). 





yO ALS ETE PALIT PT SP 


376 M. . Seaton, The chemical composition of the interstellar gas_ Vol. 111 


composed principally of hydrogen and oxygen, are more probable than metallic 
grains, and this is consistent with the suggestion that the abundance of metals 
in the gas is in substantial agreement with mean cosmic values. 


Department of Physics, 
University College, London : 
1951 June 8. 
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APPROXIMATE INTEGRATION OF THE EQUATION OF TRANSFER 
JF. B. Sykes * 
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Summary 


The value of numerical integration in obtaining approximate solutions 
of an equation of transfer, and the different methods at our disposal, are 
discussed. It is shown that although the Newton-Cotes method, used by 
Kourganoff, is better than the Gauss method, used by Chandrasekhar, both are 
inferior to a new method, the double-Gauss, discovered by the author. The 
errors in the approximate values of the source-function and the limb- 
darkening in all three methods are tabulated for various approximations, and 
illustrated by graphs. 





Numerical integration is a device to which we must turn when a function 
has to be integrated but its indefinite integral is unknown. This may arise either 
because the function cannot be integrated in elementary terms, or because it 
is itself unknown. Numerical integration consists in replacing the integral by 
the sum of the values of the integrand at a number of discrete points of the 
interval, each suitably weighted. 

If the function is known, or its graph is available, these values are easily 
found. In radiative transfer problems, however, nothing is known of the 
function /(7,.)—the intensity of radiation at optical depth 7 and in a direction 
making an angle arc cosy with the outward normal—except that it satisfies the 
equation of transfer, for instance 


pe <Iru)—4 [Men ')dp 


for a plane-parallel grey atmosphere in radiative equilibrium and _ local 
thermodynamic equilibrium (the “standard case’’) or a plane-parallel iso- 
tropically-scattering grey atmosphere. 
However, if the integral on the right can be replaced by the sum 
n 
Fel a,I(r, 145), 

where the , are a finite number of points in the interval (—1, 1), then we obtain, 
by replacing » by each yp, in turn, a set of linear differential equations with 
constant coefficients, which are easily solved to give J(7,u;). This procedure, 
which was first suggested by Wick (9), amounts to dividing the radiation field 
into distinct beams in the directions y,, as was done for two beams by 
Schwarzschild (8) in his original approximate solution of the equation of 
transfer. The solutions are of the formt 


| L, exp (—k,7) 
4 Bn i ee 
I(r, Bi) = ; [= I + LR, 


+H+O+r| . 


* Skynner Senior Student in Astronomy, Balliol College, Oxford. 
t In this and the next two equations, m is even and ,4;,..., 4n have been renumbered p-_», 





eee, Br, With p_p=—py. 


26 
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where k, is a positive root of the equation 

2 

jo-vI tp, k, 

and the L, and Q satisfy the boundary conditions 
v—1 


om 
= tak +Q=p,(t=I,...,%), 


expressing that there is no radiation incident on the surface of the atmosphere. 
These results were first obtained by Chandrasekhar (ra). 

We must now consider how the points of division and the weights are to be 
chosen. If F(u) is the function 


=2, 


ha (u— py), 
j= 
where yu, are the points of division, then 
> F(u) 
PH) = ea FOD) 


is Lagrange’s polynomial of degree less than or equal tom—1 taking the same 
values as I() for »=p,;, and 


[#0 du= 3 al) 


where 


[’ Fu) y 
F i) 1-1," 

Thus if the py; are given, the weights a; can be calculated once and for all. 
There are now three chief methods of dividing the interval (—1, 1) :— 

1. Newton-Cotes method.—In this, the oldest system (6), the interval is 
divided into equal parts, there being ‘‘ open” and “closed”’ divisions according 
as the end-points +1 are omitted or included. The weights are 


. ae He 
a= | q (os =) 


This method stands alone when Pte graph of the function is given, so that the 
values corresponding to a given » can be measured; for example, the indicator 
diagram of a steam engine. 

2. Chebyshef method.—Here the weights are taken equal, and the division 
is given by the zeros of the sili pari of 





a;= 


r exp] a 4.5x* 6.7x8 

This method is useful in the case of values equally liable to error, for the 
best result is obtained by equal weighting, but has the disadvantage that for 
high values of m some weights become imaginary. 

3. Gauss method.—Both the above methods integrate exactly with n points 
of division any polynomial of degree less than nm. Gauss (2) in 1814 inquired 
how the yu, should be chosen so as to integrate exactly polynomials of the highest 
degree possible. Since there are m roots and n weights at our disposal, this 
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degree is 2n—1. Because f—¢ is zero for ».=,;, where ¢ is defined in terms of 
S(u;), we have 


flu) - $0) = FW) Sau 


if f is of degree 2n —1, and the integration is exact if and only if 


1 
| F(u)u! du =o for 1=0, 1,..., n—1. 
om] 


It can be shown that this is uniquely satisfied by F(u)=P,,(), the Legendre 
polynomial of degree n. ‘The division points are therefore the m zeros of 
P,(u) in (—1, 1), and the weights are given by 


a;= 





r (' P(x) 
P,,'(us) ¥ Bh py; - 

This last method was the one selected by Chandrasekhar (1b) in dealing 
with the transfer problem mentioned above, and in calculating numerical values 
from his results. A variant of it, due to Radau, has been used by Kopal (3) in 
the reduction of eclipse observations of the Sun’s limb; here F(u)=(1 —p*) P; (x), 
and has zeros at » = +1, which makes it more suitable for his purposes. 

Kourganoff (4b) has recently challenged the apparently obvious superiority 
of Gauss’s formula, on the following grounds :— 

(1) I(r, ») is not representable by polynomials. 

(2) The zeros of P,() are crowded towards .=1, where J is varying slowly, 
and spaced out towards .=0, where J varies rapidly. 

In proof of (1), Kourganoff insists that J(0, ) is discontinuous at » =o and so 
cannot be represented by a polynomial. This, however, does not appear to hold 
as soon as we penetrate beneath the surface, as the many approximations to 
I(r, 2) by series in yw” or in P,(u) have shown. 

(2) is true enough in the case of even-order Legendre functions, which have 
been exclusively used by Chandrasekhar. ‘The odd-order functions have a 


zero at 4=0, but this only leads to further difficulties. ‘The function g(r), 
defined by 


Teoy=h f Mrvw)de= VP +47) 


is known to take the value 1/4/3 for s=o—the Hopf-Bronstein relation. If, 
however, the division at 1 =o is given the weight ap, it can be shown that gq, is 
given as : 

I I 

V3 V(1—4/2)’ 
and this cripples the approximation from the start. 

Kourganoff proposes instead to apply to the radiative transfer problem the 
Newton-Cotes method, with points of division evenly spaced, which at any rate 
meets objection (2). He calculates the function g(r) for various values of + in 
the second, third and fourth approximations (namely those with four, six and 
eight points in the whole interval) and also the error 


(r + q)eale op (7 + exact 
(7 - exact 





> 
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and compares them with the Gaussian values in similar approximations. 
(The exact solution for g, by which we judge the error, has been found by 
Mark (5) by means of symbolic calculus, and by Kourganoff (4a) by variational 
methods.) It turns out that the Newton method is better in the even 
approximations near the surface (0<7r<1), the essential region according to 
Kourganoff, but is slightly worse in the third. The success of Newton’s method 
is in part due to the changing sign of the error, which tends to keep it near zero; 
on the other hand, this brings about a serious difficulty, in that the error in the 
fourth approximation is greater than that in the second over r=0-4 to 0-6, a 
rather wide range, while the error in the third over this range is greater than 
either! In Gauss’s method, on the contrary, every error curve is simply a 
diminished version of the last, and one can at least be sure that by going to a 
higher approximation one is getting a more accurate result. The expression 
of the error in this form, though justifiable, also tends to conceal the relative 
inaccuracy of q,.(=@Q) in the low Newton approximations. A further objection 
to Newton’s method is the following :— 

In the course of the solution of the differential equations, it is necessary to 
express the flux integral 


1 
: I(r, o)p dy 
also as a sum over discrete values »;. Chandrasekhar puts this equal to 
n 
2 aps], 
j=1 


where the a, are the same weights as before, whereas according to the strict 
theory of numerical integration it should be equated to 


. I ft #F(x) 

= 5,1, where b,= — d 

“a. YF (uj) J-1 BBs 

If, however, F is the Legendre function, we have in fact that these two 
sums are precisely the same, since 
1 pF(u) # FY) 4, 

b =“ - 4 C d, v f ‘ 
1S” Fe) aba Foo “ones 


I [ \d 
= — =0, 
PG) ).°)* 
and the same holds for the integral 
1 
| nm 1)? dp 


which occurs in the Rayleigh-scattering case, provided the order of P is greater 
than one. The Newton method, of course, does not yield this relation. 
With these considerations in mind, a new function was tried, namely 








d" 
qn W(t I for O<p<I, 


and the same function with the sign of « changed for -1<y<o. This function, 
which is simply the Legendre polynomial transformed into the half-interval, 
has zeros at neither zero nor unity, and both the zeros and the weights are 
symmetrically distributed in the four quarters of the interval (—1, 1). In odd 
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orders, 1 =} is a zero, and the weights in even orders are half the Gaussian 
weights in half the order. We are in fact using the Gauss formula twice, just as 
is often done with the simple formulae of numerical integration to secure 
greater accuracy. It will thus be seen that the calculations are considerably 
simplified, and exactness in the radiative transfer problem is enhanced. A few 
preliminary integrations of functions such as sin‘ 7x/2, tanmx/4 and e* between 
0 and I gave excellent results, and the astrophysical calculations were then begun 
with every hope of success. Following Kourganoff, the g-function and the 
limb-darkening 
os Fe 
1(o, »)=3F (u+9+"5 “2 


ent + kp 
were calculated and the errors 


A(r+q) = Al, 
T+q Bs 


were found, the latter against the exact solution by Placzek (7). No more 
computation is needed than in either of the previous methods, and in fact (as 
indicated above) sometimes less. The objection against the Newton method, 
that the higher approximations are sometimes less accurate, holds against the 
new method (which may perhaps fittingly be called the double-Gauss method), 
but in a much more restricted range of 7, there being only one tabulated value 
where the third approximation is worse than the second. The correct evaluation 
of {Judy and {Iy*du, mentioned above, is naturally retained by the double-Gauss 
method. 

The accompanying tables and graphs show the results obtained. For q, the 
second double-Gauss approximation is better than the third Newton and 
(except for r>0-7) than the third Gauss. The new third approximation is 
much better than either of their fourths, especially for large 7. For 
limb-darkening, the second double-Gauss is better than the third Newton and 
the fourth Gauss. The new third approximation is far better than any other 
so far obtained, the relative error being less than 0-02 per cent over a great part 
of the range. The root-mean-square errors in each case were also computed, and 
emphasize the superiority of double-Gauss. 

The accuracy of the various approximations for limb-darkening in the 
range 0<p<o-I was also investigated, the maximum errors being Newton 
2°33 per cent, Gauss 3°33 per cent and double-Gauss 1-38 per cent (second 
approximations). In the third approximation the new method again outdistanced 
the others. 

The achievement of the greatest possible numerical accuracy, while not 
essential in the standard case, where the exact solution is known, is much more 
important in cases such as that of Rayleigh scattering, with the equation of 
transfer 


dl , 6 1 ’ , : , ‘ , 
pC) a 1¢7,u)— S| (3—u2) [' Hes’)du’ +(3ut@=2) [ Hrsnn'tde’ |, 
dr 16 oll -1 
and in more complicated situations, where exact solution is impossible. As 


shown above, the orthogonal property of the Legendre polynomials is a vital 
factor of success in these instances. 
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1A(r+q)l 
(7+) ox 








Fic. 1.—Errors in source-function. 
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Fic. 2.—Errors in limb-darkening. 


Second double-Gauss approximation. 
- Third Gauss approximation. 
Third Newton-Cotes approximation. 
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Fic. 3.—Errors in source-function. 








Fic. 4.—Errors in limb-darkening. 


Third double-Gauss approximation. 
Fourth Gauss approximation. 
Fourth Newton-Cotes approximation. 
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TABLE I* 
Second approximations to q(r) 

q exact qoa Ine qa Ana 
0°577351 0°57735 0°57735 0°5774 ° 
0°588236 o*58191 0°58042 0°5796 —1°06 
0°595391 0°58632 0°58342 0°5819 1°47 
0°601242 0°59057 058639 05841 1°69 
0610758 0°59866 0°59217 0°5883 —1°83 


0°627919 0°61657 060576 0°5982 —1°56 
0°649550 0°64431 0°62965 0°6154 —o°62 
07663365 0°66393 0°64974 0°6295 +006 
0°673090 0°67781 0°66664 0°6410 0°56 
0°680240 068762 0°68085 0°6505 0°63 
0°685801 0°69456 0°69279 0°6583 0°68 
0°690123 0°69947 0°70285 0°6647 0°67 
0°693534 0°70294 0°71130 0°6699 0°63 
0°696276 0°70540 0°71840 0°6743 0°57 
0°698540 0°70713 0°72438 0°6778 +051 


0°*701899 0*70923 0°73363 0°6831 +0°38 


0°705130 0°71058 0°74269 0°6879 +0°25 
0°707916 o’71119 0°75039 06918 o'l2 
0709191 0°71130 0°75 363 06931 0°07 
0°709806 0°71132 0°75499 0°6937 0°04 
0°710120 0°71132 0°75557 0°6939 0°03 
0*710270 0°71132 0°75581 06940 0°02 
0°710398 0°71132 0°75595 0°6940 +0'02 
0°710446 0°71132 0°75598 0°6940 ° 

Root-mean-square error 0°79 


* In each table, A denotes the relative error per cent. 


TABLE IA 
Second approximations to I(o, 1) 


DG NC G Ane 


0°43301 0°43301 0°4330 ° 
0°54012 0°53386 0°5278 —1°16 
0°62802 0°62413 0°6175 0°62 
071123 0°7009021 0°7038 0°28 
0°79210 0°79137 07878 —0'08 
0°87156 087171 0°8702 +0'02 
0°95009 095085 0°9513 0°08 
1°02796 1°02915 1°0314 “ O12 
1°10535 1*10684 1°1108 o°13 
118238 1°18408 1°1896 O'14 
1°25912 1°26098 1'2679 +or°15 

Root-mean-square error 0°44 
Maximum error —1°38 
Value of » for maximum error 0°04 
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TasBLe II 
Third and fourth approximations to q(t) 
Third Third 
exact dpe Qxo qa Ave Axc Aa 
0°577351 9°57735 °2°57735 °°5774 A] ° ° 
0°588236 0°58378 o'58092 05809 0°74 —I1'22 
0°595391 0°58979 0°58436 0°5843 9 O'01 1°79 
0°601242 0°59542 0°58767 0°5875 0°92 2°15 
0°610758 060560 0°59395 0°5938 —0'78 —2°54 


0°627919 0°62593  0°60777. 06080 —0'27. —2°77 
0°649550 0°65181 o'62901 0°6307 +0°27 2°42 
0°663365 0°66660 0°64404 0°6477 0°33 2'O1 
0°673090 ©00°67587 0°65487 0°6603 0°26 1°70 
0°680240 0°68219 0°66283 0°6698 0°16 1°48 
0685801 0°68686 0°66879 0°6770 0°08 1°32 
0°690123  0°69051 0°67337 0°6824 +0°028 1°20 
0°693534 0°69349 0°67695 0°6866 —0'003 I'rl 
0°696276 0°69597 0°67982 0°6898 o'o19 1°03 
0°698540 0°69806 0°68217 0°6924 —0°028 —0-96 


0°701899 0°70136 «=60°68575 + 0°6959 —0'029g —0°85 


0°705130 0°70473 0°68940 06990 —o'018 —o°7I 
0°707916 0°70784 0°69300 0°7016 —0-'003 0°55 
©°709IQI 0°70929 0°69495 0°7025 +0°003 0°44 
0°709806 0°70997  0°69604 0°7033 0°004 0°37 
©°710120 0*71029 0°69664 0°7036 0004 0°32 
0°710270 0°71044 0°69698 0°7037 0°003 0-28 
0°710398 0°71054 0°69728 0°7038 +0°002 —0'23 
0°710446 0°71057 0°69742 0°7039 ° ° 
Root-mean-square error 0°36 1°40 


TasBLe ITA 
Third and fourth approximations to I(o, 1) 
Third Third 
a4 :exact DG NC G Ane Axo Ag 


0°43301 0°43301 0°4330 ° ° ° 
©0°54012 0°53833 0°5280 —0°33 —2°24 
0°62802 0°62740 0°6149 o'10 2°09 1°85 
O°7112% O'71112 0°6979 —0'02 1°87 
0°79210 0°79218 0°7789 +O°OIo! 1:67 1°29 
0°87156 087170 0°8584 o'o161 1°51 1°10 
0°95009 0°95027 0°9371 0-0189 1°37 0°96 
1'02796 1°02816 1-o150 00184 1°26 
1°10535 I°I0555 1°0925 00172 1°16 
118238 1°18257 1°1697 0°0169 1°07 0°68 
I'25912 1°25930 1°2465 +o°oI5r —1°00 —0°62 


Root-mean-square error orl! 1°58 1°28 
Maximum error —o'69 —2'24 —2°'16 
Value of « for maximum error 0°02 ovr! 
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AN EYEPIECE INTERFEROMETER 


W. S. Finsen 
(Received 1951 April 5) 


Summary 


An eyepiece interferometer is described, intended primarily for the 
observation of close double stars. Some of the more significant features of 
the instrument are :— 


1. The interferometer unit is complete in itself, and only requires a collimator 
unit to adapt it to any telescope. 

2. It is of convenient size, with an over-all length of approximately 5 inches. 

3. As the slits are located at a pupil in a collimated beam the instrument is not 
sensitive to reasonable displacements of the star from the axis, nor does 
the slit constant depend critically on the location of the slit plane. 

. It may be ‘‘twirled”’ very rapidly between finger and thumb, so that a 
scanning procedure may be employed, of great value in eliminating spurious 
minima due to poor seeing. 

. The slit separation may be varied with the touch of a finger. 

5. The star may be examined with full aperture at will. 





1. Introduction.—In 1933 I improvised an experimental double-star inter- 
ferometer of the Anderson type (1) for the measurement of a short list of close 
pairs with the 264-inch refractor of the Union Observatory (2). Though the 
results so obtained were more than encouraging (3), the tedious operation of 
interchanging it with the filar micrometer at frequent intervals prevented its 
regular use. The obvious requirement was a combined instrument so that either 
micrometer or interferometer could be used at will with the minimum delay. 
An additional disadvantage of this interferometer was the cumbersome procedure 
for changing the slits and the necessity, in the absence of cross-threads to define 
the centre of the field, for continual checking of the projection of the slits on the 
object glass. 

In 1948 an aluminium extension tube was fitted to the micrometer, directed 
towards the object glass and provided with a dowelled stage for the reception of 
slit diaphragms. A convenient hand-hole in the telescope tube gave access to 
the stage. As an afterthought, provision was made for varying the slit separation 
by means of a Bowden cable operated from the eye end, the construction of the 
micrometer precluding the use of a more positive method of control. 

Though the interferometer could now be put into use at a moment’s notice, 
there was still room for improvement. The Bowden cable limits rotation in 
position angle to one full revolution, and the micrometer itself is not suitable for 
the rapid rotation or “twirling” so desirable in the interferometric observation 
of double stars. ‘The problem had, in fact, been narrowed down to the design 
of an external attachment to the micrometer capable of free and rapid rotation 
independently of it. 

The solution becomes obvious on realization of the fact that the interferometer 
slits may be placed at any point in the optical train. The ideal location in front 





388 W. S. Finsen Vol. 111 


of the object glass is ruled out for practical reasons. In the conventional 
Anderson interferometer the slits intersect the converging cone of rays at some 
point between the object glass and its focus. ‘This arrangement is very sensitive 
to errors in the assumed location of the slit plane relative to the focal plane, and to 
displacements of the star from the collimation axis. Finally, there is no theoretical 
objection to placing the slits in the beam of parallel rays between the eyepiece 
and the eye itself. This is, in fact, the analogue of the ideal position in front of 
the object glass and possesses comparable advantages. In practice, however, 
serious difficulties would arise in registering the slits in the minute exit pupil 
resulting from the necessarily high magnifying power. But the same effect is 
easily achieved by placing the slits at the larger exit pupil yielded by a low-power 
eyepiece, using a small secondary telescope to give the desired final magnification. 

As the slits are located at a pupil, the instrument will be insensitive to reasonable 
axial displacement of the image. It is assumed that collimation is obtained with 
the use of a converging lens; with a diverging lens (Galilean system) there is no 
real exit pupil at which to place the slits, and this advantage is lost. 

A low-power first eyepiece or collimator lens yields a large exit pupil, which 
is of advantage for minimizing the effect of errors of construction or register of 
the slits. But the exit pupil or eye-point is correspondingly more distant from the 
focal plane of the principal object glass; furthermore, a correspondingly higher 
power will be required of the secondary telescope. It is necessary, therefore, to 
effect a compromise in order to keep the dimensions of the instrument within 
reasonable limits without demanding impracticably small tolerances in the 
construction of the slit mechanism. 

In 1950 September I completed an eyepiece interferometer of this type 
(Plate 5). Since then it has been in regular use, performing under practical 
conditions with entire satisfaction. 

It is convenient to consider the whole instrument as consisting of two 
distinct units: a collimator and an interferometer. 

2. The collimator.—This is shown at A in Plate 6. The collimator lens is 
an achromatic doublet of 37} mm. focal length, mounted in a rigid housing which 
may be screwed into the eyepiece adapter socket of the Repsold micrometer. 
The flint component faces the incident light, and provision is made for accurate 
focusing, centring and squaring-on. It yields an exit pupil 2-31 mm. in diameter 
(corresponding to a power of 291) situated 4-1mm. outside the housing. 
A sleeve D rotates freely but without shake about the axis defined by the optical 
centre of the collimator lens and the centre of the exit pupil. This collimation 
must be carried out with high accuracy, but the procedure is quite simple. The 
properly squared-on lens is first focused by exarnining the image of the micro- 
meter threads with a telescope focused for infinity. It is then centred with respect 
to the axis of rotation by clamping the rotating sleeve in a vee-block, the same 
telescope being used to examine the image of an object placed at the anterior 
focus of the rotated collimator lens. The collimator housing is then replaced on 
the micrometer, and the exit pupil centred on the axis of rotation by means of 
three levelling screws not shown inthe photograph. Even witha stable instrument 
this adjustment should be checked from time to time; a convenient way of 
doing this will be described later. 

A ring B, carrying an index and reader for position-angle determination, 
slips over the rotating sleeve and is clamped to the non-rotating base of the housing. 
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3. The interferometer.—The interferometer C consists of a miniature telescope, 
focused for infinity, with an achromatic object glass of 35mm. focal length. 
Three pre-focused eyepieces are available, with effective focal lengths of 16-4, 
13-0 and 8-7mm., yielding final powers of 620, 770 and 1170. All optical 
components, including the collimator lens, have been “ bloomed”. 

The slit mechanism is located immediately in front of the miniature object 
glass. ‘T'wo slit plates E slide freely but without play in a pair of grooves. 
Springs F press a stud on each plate against opposite sides of a wedge attached 
to the plate G which moves in a slide at right angles to the direction of movement 
of the slit plates. A stud on the wedge plate engages in turn in an eccentric groove 
cut in the face of a graduated dial H which rotates about the axis of the instrument. 
This gives a large reduction factor, a half-revolution of the dial corresponding to 
a linear displacement of about 1-14 mm. in each slit plate. 


te ae —~ 
































Fic. 1.—Detail of slit mechanism. 


The slit apertures are cut very close to the opposing edges of the slit plates. 
The variable gap between the plates is optically closed by a second masking 
wedge formed by two triangular apertures in the register wedge mentioned 
earlier. At one extreme position of the dial the slit plates are withdrawn entirely 
from the optical train and a circular aperture in the register wedge permits the 
examination of a star with full aperture. 
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A tubular extension I of the interferometer carries a position circle and is a 
snug fit on the rotating sleeve of the collimator, where it is held in place by a 
spring catch J engaging in a groove. 

It goes without saying that the slit mechanism must be very accurately con- 
structed. ‘Tolerances not exceeding 0-oor inch are desirable. It is equally 
obvious that the collimator lens and its exit pupil should be centred to the same 
order of accuracy. Strictly speaking, the same high degree of accuracy is perhaps 
not necessary in the case of the remaining optical components, but I have neverthe- 
less found it well worth while, even in the case of the eyepieces, in order to prevent 
wandering of the image on rotating the instrument. ‘This ensures also that the 
centre of the apparent field will coincide with the collimation axis and therefore 
(assuming a well-adjusted principal object glass) with the point of best fringe 
visibility. 

4. The slits.—Two interchangeable sets of slit plates are available. In one 
pair the apertures are circular, and of the same size as the collimator exit pupil. 
The resultant “slits” are therefore lune-shaped, i.e. formed of two equal 
segments joined by their common chord. As the plates move inward their 
width increases with a consequent decrease of inter-centroid distance and hence 
of resolving power. The result is maximum resolving power for given light 
transmission. ‘This has proved to be a very convenient arrangement—when the 
telescope is directed to a particular star and the slits are slowly opened till there 
is sufficient light to make the fringes clearly visible, one has the assurance of the 
maximum resolving power practicable in the circumstances. 

The second pair of slit plates has rectangular apertures yielding “slits” of 
segmental shape. 

It is somewhat difficult to assess the relative merits of lune-shaped and seg- 
mental slits. ‘The lune may be regarded as composed of an infinite number of 
elementary rectangles of varying width but all centred on the axis of symmetry 
of the lune itself. In other words, the inter-centroid distance of corresponding 
pairs of elementary rectangles is constant. In the case of rectangular slits it is 
known that the slit constant depends on the inter-centroid distance alone, providing 
attention is paid to the visibility of the fringe pattern as a whole (4). This is 
therefore true for lune-shaped slits also. 

Segmental slits may also be regarded as composed of elementary rectangles, 
but in this case both width and inter-centroid distance vary. The rigorous 
computation of the fringe spacing presents an intractable problem. When I 
first used slits of this shape in 1933 I assumed simply that they were equivalent 
to rectangular slits with the same inter-centroid distance. In this I have been 
followed by Wilson (5), who further points out that in the limiting case (never 
realized in practice), when the segments are semi-circles, this procedure gives a 
“slit constant’’ only 3-4 per cent smaller than the theoretical resolving power of 
an unobscured circular object glass. For infinitely narrow segments the error 
is, of course, zero. In computing the slit constants for the segmental slits I 
have made an adjustment for this error based on linear interpolation. 

Apart from the slight uncertainty in the value of the slit constant, segmental 
slits seem to have some theoretical advantage over lune-shaped slits. For the 
same inter-centroid distance they transmit about 6 per cent more light, and the 
diffraction pattern is clearly narrower and probably longer. This was confirmed 
in a laboratory test. The shorter fringes were apparently brighter and easier to 
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see, and there was noticeably less contrast as regards both shape and brightness 
between central and extreme fringes. As a result both minima and “zig-zag 
effects’’ (2) were decidedly easier to recognize. A similar test was made under 
practical conditions at the telescope ; variable seeing prevented a clear-cut decision 
but the conclusion was, on the whole, in favour of segmental slits. 

With the present instrument used in conjunction with the 26}-inch refractor 
the slit constant can be varied between 0”:087 and 0”-132 with lune slits and 
between 0”-087 and 0”-143 with segmental slits. In both cases the smaller values 
correspond to infinitely narrow slits and are not realizable in practice. 

A mask with a circular aperture of the same size as the exit pupil formed by 
the collimator can be placed immediately in front of the slit apertures. This is of 
doubtful value with a well-adjusted instrument and is usually dispensed with. 
It has the disadvantage of stopping-down the instrument as regards visibility of 
the micrometer threads. 

5. Adjustment of the interferometer.—It is assumed that the optical components 
have been centred relative to the rotation axis as described earlier. The collimator 
exit pupil must now be centred, at the same time checking its registration with 
the slits. As this operation should be repeated from time to time, a special 
attachment has been provided for the purpose. This consists of an extension 
tube fitting over the eyepiece, with the help of which one of the remaining eyepieces 
can be focused on the final exit pupil. The slits are thus seen projected on the 
principal object glass, making it easy to centre the collimator exit pupil by means 
of the levelling screws mentioned earlier, at the same time checking the centring 
of the slits. (This simple procedure is useful also for checking the slit projection 
of Anderson interferometers.) 

The position-circle zero is easily determined by directing the telescope at a 
bright star with narrow slits and aligning the elongated diffraction pattern with a 
micrometer thread. 

With a well-adjusted instrument the fringes should have maximum visibility 
in all position angles at the centre of the field. Near the edge of the field this will 
in general no longer be true. When the fringes are radial their visibility will be 
but little impaired, but on rotating through 90° they will become more or less 
invisible. The effect is indistinguishable from genuine duplicity, but this 
erroneous conclusion is easily guarded against by keeping the star in the centre 
of the field. 

It may happen, however, that the position of maximum fringe visibility does 
not coincide with the centre of the field. The effect is noticed with both the 
Anderson interferometer and the present instrument. On the assumption that 
the star is zenithal (for atmospheric dispersion has the same result) and that faulty 
adjustment of the interferometer has been eliminated as a possible cause, the most 
probable explanation is relative decentring of the object glass components, re- 
sulting in a prismatic effect. The fault is not cured by a reasonable squaring-on 
adjustment of the object glass, nor by moving the eyepiece laterally; the fringes 
remain invisible unless the direction of vision is more or less oblique—the pencil 
of rays passing through the eyepiece is thereby decentred so that a compensating 
prismatic component is introduced. If the star is placed centrally in the field it is 
easy to determine the sense in which the components of the object glass are 
decentred by placing a prism of small refracting angle between the eyepiece and 
the eye and noting the orientation that gives improved fringe visibility. The 
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object glass clearly requires prismatic correction in the sense indicated by the 
prism, making proper allowance for any ray inversion. 

It may be remarked in passing that this procedure provides a very simple 
and sensitive method of adjusting the relative centring of object glass components. 
For this particular purpose the slits may be cut in cardboard and placed in the 
ideal position in front of the object glass, thus eliminating other extraneous 
factors. 

It is also of interest to note that it is fallacious at large zenith distances to 
compare the fringe visibility of a suspected double with that of a neighbouring 
star as control, unless the precaution is taken of examining both stars at exactly 
the same point in the field. Fringes invisible at the centre of the field as the 
result of atmospheric dispersion may be readily visible nearer the edge with 
oblique vision. But the deliberate use of this fact to compensate for atmospheric 
dispersion cannot be recommended; it is too uncertain and it introduces the 
possibility of other off-axis errors (especially in the case of the Anderson inter- 
ferometer). ‘The safest procedure is the introduction of variable dispersion 
between the eyepiece and the eye by means of a rotary or Risley prism (6). 


Union Observatory, 
Johannesburg : 
1951 March 27. 
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PHOTOELECTRIC RECORDING OF STELLAR OCCULTATIONS 


H. E. Butler 
(Received 1951 February 26) 


Summary 


A photo-cell is mounted near the focus of a 15-inch reflecting telescope 
and the photo-current, corresponding to the intensity of starlight falling 
on it, is displayed as the vertical movement of one of the spots of a 
double beam oscilloscope. This deflection is recorded photographically on 
horizontally moving paper. Seconds marks obtained photoelectrically 
from the Observatory’s Mean Time pendulum clock are shown on the 
second spot of the cathode ray tube. When a star is occulted at the dark 
limb of the Moon the change in deflection is clearly visible and its time is 
easily measured with respect to the clock marks to about ;4; of a second. 

The same method is used to determine the clock errors, the B.B.C. time 
signals being displayed directly on the tube instead of the star. As a result, 
the clock error with respect to the time signals can also be determined to 
about ;4, of a second. 

If several measures of clock error are made before and after the 
occultation, it is believed that the time of the occultation can be determined 
to within 1 or 2 hundredths of a second of the time defined by the time signals. 





1. Introduction.—The visual method of timing the occultations of stars by 


the Moon gives a time which is not only late by the personal error of the observer, 
but which is, in addition, only accurate to about + of a second. The observer, 
using a telescope of reasonable aperture, starts a stopwatch when he sees a star 
disappear (or when he sees it appear, whichever the case may be) and then compares 
his watch with a carefully rated clock. 

Referring for a moment only to disappearances at the dark limb, which are the 
occultations most frequently observed, errors can arise in several ways: 

(a) The observer will know the expected time of the occultation with an 
accuracy which will depend upon the geographical position of his observatory 
with respect to those observato=es for which occultations are given in the Almanacs. 
This uncertainty may amount to as much as one or two minutes. The majority of 
observable occultations occur when the Moon is between 7 and 14 days old and 
when, as a result, the dark limb will be invisible. In most cases, therefore, the 
observer will have no warning of the occultation and will start to observe the star 
well before it is actually occulted, and his eye will be appreciably fatigued when 
the star eventually disappears. The suddenness of this disappearance is always 
surprising, and frequently one knows that up to half a second has elapsed between 
seeing the star disappear and starting the stopwatch. 

In comparing the watch and the clock, since the seconds hand of the clock 
can be seen and the clock-ticks can be heard, the greater part of the personal error 
is avoided. Nevertheless, random errors of up to about 4, of a second are still 
to be expected. 

(6) The stopwatch itself is seldom a precision instrument. It may even have 
a small error over the few minutes that it will be running, but the biggest error 
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will usually be that it moves in jerks of + of a second at a time, which, coupled 
with the difficulty of reading, introduces an error of up to } of a second. 

(c) The master clock will always have an error, and we must consider how 
accurately this error can be measured. Our own procedure is either to look 
directly at the clock and compare it with the B.B.C. time signals heard from a 
radio set in the same room, or to do the same thing using a stopwatch as an inter- 
mediary. In either case the error will be up to + of a second. These errors 
can be reduced to some extent by making several readings before and after the 
occultation and graphing the result, but even so the final error will be of the 
order of +, of a second. 

(d) If a chronograph is in use, (6) is eliminated and both (a) and (c) are 
unaltered. We introduce, however, the troubles that can arise when contacts 
are in use on an accurate clock and also introduce electrical and mechanical lags, 
etc., in the pens. 

When observing appearances at the dark limb, the errors due to (d), (c) and 
(d) are unaltered, but that due to (a) is considerably increased. 

2. Method.—It appears that the following method enables one to time an 
occultation with respect to the broadcast time signals with an accuracy of I or 2 
hundredths of a second. 

A 1P21 photomultiplier tube has been mounted about 2 inches behind the 
focus of the 15-inch reflector. ‘There are changeable diaphragms at the focus. 
The output from the photo-cell is fed through an amplifier to give vertical 
displacements to one of the spots of the double-beam cathode-ray tube of an 
oscilloscope. A camera records these deflections on horizontally moving 


photographic paper, the spots themselves moving only vertically. 

The second spot of the cathode-ray tube is given a deflection once a second by 
the Observatory’s Mean Time pendulum clock. As some trouble had been 
experienced with the electrical contacts upsetting the rate of the clock, it was 
decided to use a photo-cell and a light beam, and so to impress seconds marks 
on the oscilloscope record without interfering in any way with the clock. 

The arrangement is shown in Fig. 1. 


Cc 


Fic. 1.—Arrangement for taking off seconds marks from clock. 
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Fic, 2.—Double beam oscillograms. 


Lower trace in each case shows seconds marks from Observatory clock. 
Upper trace shows: 
(a) and (c): Broadcast six pips at 19!" oo” and 22" oo” respectively. 
(6): Trace of star y Cancri, with occultation at point A. 


H. E. Butler, Photoelectric recording of stellar occultations. 
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Consider the pendulum in its vertical position. A small plane mirror, 
2x6mm. (M,) (actually a galvanometer mirror) was fixed to the point of the 
pendulum bob, with its longer side horizontal, and set at an angle of 45° to the 
horizontal. A small lamp L, with lens L,, is mounted outside the glass clock-case 
and directs a beam of light horizontally so that the mirror M, is illuminated. 
Beneath M, is fixed a diaphragm with a slot S which is vertically under M, and 
the same size as the projected area of M,. A portion of the original beam is 
reflected downwards by M,, through the slot S and on to another mirror M,. 
This mirror directs the light back horizontally through another lens L, and out 
through the case to fall on the cathode of a photomultiplier cell C. 

As the pendulum swings, the cell is illuminated only as M, passes over S. 
The output current of the cell is fed through a cable to the telescope dome where 
it is displayed as the vertical deflection of the second spot of the oscilloscope. 
Each swing of the pendulum thus gives a short time-mark similar in shape to that 
of achronograph. The duration of the mark is about 0-05 second. The lower 
traces of each of the three examples in Fig. 2 (Plate 7) show these marks when 
recorded with the relatively slow paper speed of 0-4 inch per second. 

3. Accuracy of rating of clock.—Time signals are broadcast at most whole 
hours by easily received, medium wave-band transmitters of the B.B.C. from 
05" 00™ to 23" 00™ and other times can be filled in, if necessary, with short-wave 
transmissions. ‘The clock error is determined photographically at convenient 
intervals before and after the occultation, in the following way. Seconds are 
still shown on one spot of the oscilloscope, but the amplifier of the other spot, 
instead of being connected to the photo-cell on the telescope, is connected to the 
speaker terminals of a normal radio set. Sound is then displayed on the tube and, 
in particular, the six pips can be clearly seen and recorded. Examples are to be 
seen in the upper traces of (a) and (c) of Fig. 2 (Plate 7). 

Fig. 3 shows clock errors obtained in this way on two occasions before and two 
occasions after an occultation, together with the best that could be obtained with 
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Fic. 3.—Clock correction near time of occultation, 
Ordinates: Clock correction (seconds). 
Abscissae: Time (1951 February 19). 
The arrow indicates clock correction at time of occultation = +- 19°87 seconds, 
x : Clock correction by stopwatch. 
©: Clock correction photographically. 
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a stopwatch. Traces (a) and (c) of Fig. 2 (Plate 7) were the actual traces which 
were used to obtain the points at 19" 00™ and 22" oo™ respectively. 

It will be seen that a very smooth curve can be drawn through the photo- 
graphically determined points, the scatter being negligible. If we have a reason- 
ably good clock and make frequent comparisons with the time signals, the time 
at any moment can in this way be determined to within ;+, of a second of the time 
as defined by the time signals themselves. 

4. Example of a recorded occultation.—Fig. 2(b) shows (Plate 7) the trace of 
the star y Cancri, of magnitude 4-7, as it was occulted by the Moon (13:5 days 
old) on 1951 February 109. 

Even with a focal aperture of 0-5 mm., considerable trouble was experienced 
with glare from the Moon, which was nearly full. As the Moon approached the 
star, the deflection of the spot steadily increased as a result of this scattered light. 

The sudden decrease in deflection as the star was occulted is apparent as the 
downward movement of the trace at A. Prior to A (on the left), the scintillation 
of the star can just be seen on the trace, while after A only instrumental noise 
is present. 

The speed of the recording paper was 0-4 inch a second and, needless to 
say, no idea of stellar diameter could be obtained. The break at A is sufficiently 
sharp for its time to be measured as 0-47 + 0-01 second after the nearest seconds 
mark on the left. The occultation was also recorded by stopwatch in order to 
relate the fraction of a second so obtained to the appropriate whole minute and 
second. ‘This enabled the time of the occultation to be written as 19" 37™ 518-47 
referred to the Observatory clock. Allowance for the clock error as read from 
the graph of Fig. 3 gave the final corrected time of the occultation, as obtained 
photographically, to be 19"38™ 115-34. Using a stopwatch both for the deter- 
mination of clock error and for timing the occultation itself, gave a time of 
19" 38™ 118-4, 

5. Conclusions.—The method described should give occultation times for 
disappearances at the dark limb as accurately as the broadcast signals warrant. 
It will be used on the brighter stars in future occultation programmes, both to 
obtain accurate times and to see what errors are present in normal observations. 
Actually, in order to attempt measures of stellar diameter at the same time, higher 
recording speeds will be used, although the resultant greater timing accuracy 
of the occultations is hardly justified. 

It is hoped to try the method for reappearances at the dark limb. It will be 
necessaty to centre the star in the focal diaphragm before it disappears behind 
the bright limb, and to rely on the telescope driving accurately for the whole time 
that the star is out of sight behind the Moon. Since this time can be as much as 
one hour, this may prove impossible. 

The method is not suitable for appearances and disappearances at the bright 
limb. 


Dunsink Observatory, 
Co. Dublin : 
1951 February 23. 
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Summary 


Models consisting of a hydrogen envelope surrounding a core of higher 
molecular weight, with opacity throughout the star due to electron scattering, 
are examined. It is found that considerable (though only finite) extensions 
of radius can occur with a moderate fraction of the mass of the star in the 
envelope. The evolutionary tracks in the Hertzsprung—Russell diagram are 
examined, and good agreement with observation is found. 





1. The present paper constitutes a further step in the systematic examination 
of the red giant problem which is the purpose of this series. In Paper I (x) the 
general features of the problem were discussed, and it was shown that no 
physically plausible law of opacity could lead to a homogeneous star that generates 
energy by the thermonuclear processes and possesses the large radius characteristic 
of red giants. Paper II (2) began the examination of stars possessing a single 
discontinuity of composition. The special opacity law « ~Z(1 +X)°75p°757-*5 
was assumed to apply, and it was shown that large (and even mathematically 
infinite) radii could occur provided the change in composition was reasonably 
large and provided it took place in a certain part of the outer layers of the star. 
The model is then the Hoyle-Lyttleton model (3, 4) of a red giant, with a special 
law of opacity. 

The corresponding problem for the case of opacity due to electron scattering 
forms the subject of the present paper. It seems likely that in actual red giant 
stars the law of opacity is intermediate between the case investigated in Paper II, 
where x ~Z(1 +X)°75p°757-35, and the electron scattering case (4) in which 
«~(1+X). It may then be presumed that the structure of an inhomogeneous 
star with intermediate law of opacity is intermediate between the two cases. 
This will be especially likely if (as is the case) the investigation reveals that there 
are no really substantial differences between the two models, at least as far as the 
red giant problem is concerned. The problem of the structure of an 
inhomogeneous star with opacity due to electron scattering is, in any case, a 
subject of considerable interest. 

2. The problem (first examined by Hoyle and Lyttleton (4)) may be stated 
as follows. The sphere of radius r=r, divides the star into two parts, each of 
which is homogeneous in composition. The abundances of hydrogen, helium, 
and other elements in the interior will be denoted by X;, Y;, Z;, respectively, 
and the mean molecular weight there will be called yp, The corresponding 
quantities in the exterior region (r>r,.) will be called X,, Y,, Z, and p, 
respectively. In accordance with the view that the difference in composition 
is due to the thermonuclear conversion of hydrogen into helium we shall assume 
that Z,=Z,. Since Z is taken to be very small in accordance with modern views, 
and is of no importance for the process of electron scattering, we shall put 
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Z,=Z,=0. The material in the outer region is taken to be accreted interstellar 
matter in which the helium abundance is likely to be very small. Hence we shall 
put Y,=o and therefore X¥,=1, »,=}. The set of stellar models constructed 
therefore forms essentially a two-parameter series, the hydrogen abundance in 
the interior portion being one parameter (Y; and py, being determined by it) 
and the mass ratio M,/M, being the other one. 

The mode of energy generation in the large stars (in which alone electron 
opacity predominates) is almost certainly the Bethe cycle. This process is highly 
temperature sensitive. ‘Therefore effectively all the energy generation will take 
place in the inner part of the star (r<r,) provided that neither M,/M, nor X, is 
too small. The first possibility concerns a case of little interest, the second may 
be shown to require X, to be less than about 0-or. In the models to be considered 
the assumption will always be made that all the energy generation takes place 
in 7 <r, so that the model does not apply if X; is less than o-or. The case X;=0 
will be considered but must be understood as a purely mathematical construction 
that can be approached fairly closely but not reached in reality. 

Owing to the temperature sensitivity of the Bethe cycle the innermost part 
of the star will be in convective equilibrium, Denoting the radius of the 
convective core by 1,, it is clear that r,;<r, since otherwise rapid mixing would 
occur. It will be assumed that all the energy generation occurs for r<7;. 
This assumption greatly simplifies the treatment, but it also must be closely 
approached in any actual star. 

Radiation pressure is not likely to be a factor of great importance in the 
models to be discussed in the paper, and will be neglected. 

3. The outer part of the star (r>r,) is either wholly or partly in radiative 
equilibrium according to the values of M,/M, and y,;/p,. It can be shown that 
convective regions occur there only for rather small values of M,/M,, a case of 
little interest for the red giant problem. Just as in Paper II, the case in which 
there is a convective zone in the region r>r, will be excluded. The treatment 
is then very similar to that used in Paper II. The law of opacity in the radiative 
region in the present case is 

«=K(I+X) (1) 
where x, is a constant, which in c.g.s. units has a value of about 0-19. 
The energy transport equation is therefore 
aT _ _ 3Ko(1+X)Lp (2) 
dr 16zacr*T* * 
We define h-variables as in the previous papers in this series (1, 2), that is 


_ 40Prt Pr R. d (log p) 
S= Gu 2= Gump N=*- Tog Py (3) 
Then the hydrostatic condition is expressed by 
SdQ_ S+N-Q (4) 
QdS 1+2S-—4Q’ 4 


while in the radiative region, the energy transport equation is by (2) and (3) 


SadN_ S+1-4N ts) 
N dS ~1+25-—40° 5 


In the convective core N=0-4. 
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At the interface r=r,, the quantities P, 7, M and L are continuous, and 
hence S is continuous. Since P and T are both continuous, p/u is continuous 
and hence from (2) 

I dT 
w(i+X) dr 
is continuous. If we denote the values of the discontinuous variables just 
outside the boundary r=r, by the suffix e and those just inside r=r, by the 
suffix i we then have 
Pac _ Be, Qre _ Hi, Nee -_ 1+X, (6) 
Pose” Og Me’ SNe 14+ 

The star may be divided into three zones as follows: 

(a) r,>r>rz. An envelope in radiative equilibrium with composition 
(X,., Y., Z, 3 -) (taken to be (1, 0, 0 ; 4) respectively). 

b) r,>r>r,. An intermediate zone of composition (X,, Y;, Z ; @,) in 
radiative equilibrium. 

(c) r,>r>0. An inner convective zone having the same composition as 
zone (b). 

Zone (c) always exists and so the range of cases arising depends on the 
relative extents of zones (a) and (6). In order that the arrangement shall be 
stable, we must have »,<p,. In the case in which zone (a) vanishes, the star is 
homogeneous and the corresponding solution will be referred to as the standard 
solution. 

In zone (5) equations (4) and (5) apply and we have a second-order system. 
Only a one-parameter family of solutions of this system is of interest in zone (d) 
since at S=S, no change in composition occurs and consequently all the 
variables are continuous at S=S,. We have to consider in (d) therefore only 
those solutions of (4) and (5) that fit smoothly on to the solution for the core. 
This is the Emden solution 

O=Qz (S;0°4), N=0-4. (7) 

We therefore require to investigate the solutions of (4) and (5) for r>r,, and 

this condition implies S<S,. The value of S, may be taken as the parameter 

specifying the members of the family. A particular member of the family is the 

standard solution in which S,=S{=1-03, which runs to the surface 
S=Q=0, N=} (5). This standard solution will be written 

Q=0*(S), N=N*S). (8) 

As in the problem solved in Paper II, the other members of the family lie close 

to the standard solution, provided that the value of N does not differ by a large 

amount from the standard value. These solutions diverge very rapidly from the 


standard once the N values deviate markedly. A higher approximation to the 
solutions can be obtained as follows. Consider 





@ 1+5S—40 
obtained from(4)and(5). It will be seen that the rate of change with S of log( NQ*) 
depends on Q and S only. It follows that if at any point S=.S, of the zone (5) 
solution, N—N* is of the same order of magnitude as N*, then at that point 
and for all S>S,, Q is close to Q*. By virtue of the first approximation the 
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right-hand side of (9) is very close to its standard values and hence a higher 
approximation is given by the relation 


NQ*= N#Q* (10) 
provided that N—N* does not differ by a large factor from N*. Equation (10) 
holds within a few per cent if this condition is satisfied. No greater accuracy 
is likely to be required in zone (6). 

In zone (a) the solutions also form a one-parameter family of solutions of (4) 
and (5) in the (.S, Q) plane, since they must all run to the surface S=Q=0, N=}. 
Near the surface, these solutions are given approximately by 

O = AS"* + O(S5*), N=}+ O(S) (11) 
where A is a constant. A may be used to label the members of the family, which 
has precisely one member in common with the family (b), namely the standard 
solution with A= A*. Asin the problem discussed in Paper II, we are interested 
only in solutions with A > A*. 

The similarity with the problem in Paper II breaks down when we consider 
the curve described there as the spiral F. In the present case, the spiral point 
is the point S=o, Q=N=}. The spiral F, if it existed, would also start from 
this point. In fact no such solutions exist, and F may be considered to reduce 
to the point S=o, O=}, N=}. 

All members of the family with A greater than some value A (A <A*) run to 
the point S=o, O=N=} (cf. Paper I). Those with A<A run off to infinity. 
The curves fill the entire positive quadrant of the (.S, Q) plane and there is one 
and only one member of the family through each point of this quadrant. 

A definite value of N therefore corresponds to each point of the (S, Q) plane, 
namely the value taken by N at that point by the unique member of the family 
through the point. Curves of constant N can hence be drawn in the (S, Q) 
plane. Reference to Fig. 1* shows that the variation of N is slow and smooth. 

As a solution is followed into the star from S=Q=0, N=}, the expression 
dS/(1+2S—4Q) is always positive and so are N and S. Therefore, N will 
increase or decrease according as S+1—4N is positive or negative. It follows 
that N cannot first reach any given value N,>} unless S>4N,—1. Hence 
convective instability cannot first occur for S<o-6. It will be shown later that 
large radii occur only for values of S, much smaller than 0-1, and hence the 
case we are discussing now is the one of interest. 

Just as a definite value of N can be associated with each point of the plane, 
so also can a definite value of M/M, and r/r, be associated with each such point. 
These quantities are given by the quadratures 

log M, = { nie leah (12) 
M = Jot+2S—4Q” 
OdS Be ee 
=; * ee | ee (13) 





The definite values are obtained by integration along the solution through 
the point. Thus curves of constant M/M, and constant 7/r, can be drawn on 
the diagram representing the envelope solutions (Figs. 1, 2). The curves of 
constant M/M, start on the S axis with M/M, equal to (1+2S)¥*, They all 
run to the point S=o, Q=}. The curves of constant 7/r, start on the Q axis 
with r/r,=}—Q and run to infinity in the (.S,Q) plane. The locus r/r,=o in 

* The figures given in this paper summarize the results of eight integrations, 
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this problem consists only of the line S=o for Q >}, and there is no curve 
corresponding to the spiral F. It follows that no stars of infinite radius can 
be constructed with this law of opacity since the equations (6) make the points 
on S=o of no interest as points of fit. 

In constructing a complete stellar model, the solution for zone (a) has to 
be fitted on to that for zone (b). At the fitting point, S, «Q and N/(1+X) are 
continuous. The composition has been discussed in Section 2. As stated there, 
we shall take Z,=Z,=o and X,=1, Y,=0. It follows that the maximum 
discontinuity possible occurs when Y,=1 and X,=o0. This condition would 
be produced when all the hydrogen in the region r<r, had been converted into 
helium in the energy production processes. 








Ns0:275 ‘phees 
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Fic. 1. 
Thick lines : Solutions (A=0-218, 0°250, 0°270, 0°338, 0°42)). 
Thin lines : Curves of constant N. 
Dashed lines : Curves of constant r/rs. 
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Fic. 1 (a).—Detailed diagram of neighbourhood 
of S=o0, Q=0'25 (solutions drawn for 
A=0°338, 0°429, 0°483, 0°593, 0°695). 
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4. In the envelope, the solution will be given by one of the curves (a) already 
described. At S=S,, both N and Q jump in accordance with equations (6) to 
one of the curves (b). For any S;,, the curves (5) relate Q,, to N2,. Therefore 
for any given change of composition, there corresponds to each value of S, a 
single set of values (Q.,, N2,), (Q2;, Nz), such that the points (S,, Q,,, N2,) and 
(S2, Qo No) lie respectively on a curve (a) and a curve (6), and such that 
equations (6) are satisfied. ‘The locus (S,, Q,,) followed as S, varies will be 
termed the “fitting curve”’. The point at which the fitting curve intersects 
any curve M,/M,=constant, gives the correct value of S, and Q,, for a model 
having a given change in composition and a given proportion (M,— M,)/M, 
of the mass in the envelope. 

To construct the fitting curve accurately would involve performing a 
considerable number of integrations both of solutions (a) and of solutions (8). 
However, the fitting curve can be constructed approximately with very little 
calculation by making use of the fact that the zone (4) solution lies (as far as Q is 
concerned) close to the standard solution. We consider the ratio Q,,/Q,.° and 
assume (10) to hold exactly. Putting X,=1 as usual, we have 


Qoe _ Dae Dai _, x) a N$\"* (N2.\"4 
of 0, 08 7*\N,) ~7\N) LN, 


2 \"4eNs\i4 
~n(rex) (x) a) 


Thus the fitting curve can be obtained by scaling up the standard solution 
by the factor on the right. This factor involves the ratio Nz'/N.,, which is 
determined completely only when the fitting curve is known, but is in fact 
always close to unity. The following limits can be placed on N.*/N2,. The 
curves of constant N are very steep for small S, and as S tends to zero, N2’/No, 
approaches unity. It is found from computation that, provided no convective 
zone occurs for r>r,, the curves of constant N are always sufficiently steep and 
sufficiently widely spaced for N.*/N,, to exceed 4(1+X,). Thus we have the 
inequality 

21, < Or.) OF <2ufh(1 +X)" (15) 


Q,,/O% therefore lying within fairly close bounds. If we take for a given S, a 
Q,, satisfying (15) we can then find the value of N."/N.,. Using this we can find 
a second approximation to Q,,/O0%. This can be continued until a sufficiently 
accurate approximation to the fitting curve is obtained. The second 
approximation is generally sufficient. 

If the fitting curve is drawn on the (S,Q) diagram for the envelopes 
displaying the curves M,/M,=constant (Fig. 2), then, as in Paper II, we can select 
the curve corresponding to a stellar model having an assigned fraction of the 
mass in the envelope. The point (.S,,Q,,) is the point at which the fitting curve 
intersects the curve M,/M,=constant, and the values of N,, and r,/r, may also 
be read off (Fig. 1). The quantities Q,, and N,, can be calculated from (6) and 
the intermediate solution (+) constructed if the standard solution is not 
considered to be sufficiently close to the actual solution. 

5. It is of great importance to have some measure of the radius of a model 
so that one can compare the radius with that of the standard model. Such a 
measure must contain the information given by the Bethe law that the central 
temperature of all models considered must be approximately the same. 
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Consider the quantity 


BT, r\_ °  NdS 
tog (“Tt om) = 982+ |) sepas apy’ (76) 





the point (S, Q) being chosen to be a point of the solution where r<r,. Since 
for r<r_g we have Q and N approximately equal to their standard values, it 
follows that log Q and the integral on the right-hand side are close to their 
standard values. It follows that for stars of the same central temperature and 
central composition, the value of 7,/M, is close to its standard value. 











0-5 


Thick lines : Fitting curves (u,/p,=1, 4, 4, 2, 5, $). 
Thin lines : Curves of constant M/M,,. 
Dashed lines : Curves of constant E. 











+S 
Fic. 2 (a).—Detailed diagram of the neighbourhood of Smo, Q=0'25. 


We now introduce a special form of (16), namely the invariant D used in 
both Papers I and II. 
', RT, 


D= Gia (17) 





404 C. M. Bondi and H. Bondi Vol. 111 


It again forms a convenient quantity to measure the radii of various models. 
In h-variables, we have 
_ (%» (Q—S)dS ° NdS 

log D= | Sasa) * 8%* |e satay | 
where (S,, Q,) is an arbitrary point of the particular solution. If, as is 
convenient, we choose (S,, O,)=(S., Q.,), we know that the second and third 
terms in the expression on the right-hand side are almost the same for a 
non-homogeneous model as for the standard model. If D* is the value of the 
invariant (17) for the standard model, it follows that the ratio D/D* depends 
only on the changes in the first integral in the non-homogeneous model as 
compared with the standard model. This ratio is clearly of great importance, 
and is as before called the extension factor and denoted by E. It is the ratio of 
the radius of any model to that of a homogeneous star of the same mass, central 
temperature and central composition. If the models we wish to compare do not 
have the same central composition, allowance can be made for this by multiplying 
E by an appropriate factor (the ratio of the central molecular weights). 

We have 


log E=| 








‘S% (O-—S)dS Ss (Q*-—S)dS 
0 S(Ii+2S—4Q0) Jo S(1+2S—4Q*) 
< .: (O—Q*)(1-—2S)dS 
=), Sat25-40\a+2S—40%) 79 
We see, as in Paper II, that from consideration of the signs of the factors in the last 
integral, no larger values of E can occur to the right of S =o-5 than occur to its left. 
In practice, we calculate E most easily by comparing the ratio Mr,/M,r at 
the point (S,,Q,,) with its value for the standard solution at S=S,. In this 
way we compare the radius of a non-homogeneous star with the radius of the 
standard-solution star having the same mass, central temperature and central 
composition. There is one value of E corresponding to each point (S;, Q2,) 
so that curves of constant E can be drawn in the envelope diagram (Fig. 2). On 
the standard solution E=1, while E is infinite only at points on the Q-axis for 
Q>}. The curves of constant E start on the Q-axis with E=(1-—4Q)"". Fora 
definite value of the ratio y,/u, and a definite value of M,/M, the corresponding 
E value may be read off the diagram containing the fitting curve and the curves. 
of constant M,/M,. ‘The relevant E value is that corresponding to the point 
at which the curve representing the ratio M,/M, intersects the fitting curve. 
Since the only points at which E£ is infinite lie on the Q-axis above 0 =}, it 
follows that for no finite value of ,/4, can an infinite extension factor be obtained. 
For, finite values of y,/u, give rise to fitting curves which reach the Q-axis only 
at the origin. If we restrict ourselves to values of .,/u, not greater than 8, that is,. 
the case of pure hydrogen for r>r, and pure helium for r<r,, then the maximum 
value of E obtained appears to be in the neighbourhood of 36. If for given 
14/h,. we plot E against M,/M,, we find that as M,/M, decreases, that is, as the 
mass in the envelope increases, E first increases, then reaches a maximum and 
then decreases again as further mass is added to the envelope (Fig. 3). 
6. We now have to deal with the values of the luminosity for these: 
non-homogeneous models. With the opacity law used here 
a (3) oa 3Ko(t+X)L (20) 
S G/ (47)? 4acM* 











at al] points. 
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Applying this at the surface, we have the mass-luminosity relation 
R\* 1 3K,(1+X)L 

ae ve, nine od 2 
sels (73) (4m)? 4acM* (21) 
With this opacity law, we see that the radius does not enter into the relation, 


so that the relation between mass and luminosity does not depend, as in the 
case discussed in Paper II, on the central temperature. 

















Beyond this 
curve 

|__ convection 
sets in 











Fic. 3.—Variation of extension factor E against M,/M, for various ratios y,/m,. 
Vertical scale: 1—E-*!*, 


By computation, it is found that a semi-empirical relation connects the 
luminosity of a non-homogeneous star with that of a certain homogeneous star. 
Equation (21) shows that assuming the same opacity law to apply to both stars 
and that their masses are equal, the luminosities are proportional to A*y,4(1 + X,)"}. 
If we take a homogeneous star with the same mean value of the composition as 
the non-homogeneous star, then it is found by computation that A4y,4(1 + X,)-} 
is almost exactly the same in the two models. This property appears to hold 
for all non-homogeneous models with this opacity law, whether the extension 
factors are large or not. 

7. In order to obtain a clearer view of the evolutionary significance of the 
theory discussed in this paper, the application of the theory to actual stars will 
now be considered. 

It will be seen from Fig. 3 that appreciable values of the extension factor 
only occur if 0-95>M,/M,>0-4. This implies that the discontinuity in 
composition occurs not at the surface of the convective core but somewhere 
in the outer layers of the star. This is a characteristic feature of the Hoyle- 
Lyttleton model of red giants. It can arise only if a well-mixed star that has 
converted a considerable fraction of its hydrogen into helium suddenly begins 
to accrete interstellar hydrogen at such a rate that it cannot be mixed in with 
the rest of the star. The required conditions are therefore that the speed of 
accretion should exceed the speed of mixing, which in turn must exceed the 
rate of helium formation (6). The currents recently discussed by Sweet (7) are 
presumably responsible for mixing the helium produced in the core throughout 
the original envelope of the star. 
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Assuming these conditions to be satisfied, the evolutionary sequence can now 
be traced.* A star of mass M, radius 7 and luminosity L, consisting of almost 
pure hydrogen, evolves through the process of helium formation. The rate of 
mixing is taken to be high enough to ensure a homogeneous composition. Then 
the new values of M, r and L are related to M, r, L by 


ait EC ey pan ets 
M=M, L= cx roma yt- (22) 


In order to consider the variation of the radius we have to consider the 
process of energy generation in some detail. For stars sufficiently massive to 
make electron opacity predominant the process of energy generation will be 
highly temperature-sensitive (Bethe cycle). The assumption of a constant 
central temperature will therefore be a fair approximation, and will be adopted 
since it greatly simplifies the analysis. Then since the invariant 

a 2 
wy, GM, 
has a fixed value for a homogeneous star, 


r= 2ur. (23) 
At some stage this slow development ceases to control the speed of evolution, 
and accretion is supposed now to become the most important factor. Let the 
fraction of the original mass added by accretion be v so that in the subsequent 
development 
M=(1+v)M. (24) 
Since M is the mass of the part of the star having molecular weight p, the ratio 
M,/M,=1/(1 +»). 
Since », and X, are taken to be constant throughout the subsequent develop- 
ment, the luminosity will vary according to 
L~A‘M? ~A(1 + v)?. (25) 
But for very small v the luminosity must be unchanged by the thin hydrogen 
envelope. Hence, by (22) and (14) (since, for very small S,, 4/A*=Q,,/Q}) 


; = (js) ( +r). (26) 


Also, for constant 7, and y,; 
r~EM~E(1 +1). 


Combining this with (22) we have by continuity at v=o 


© = 2y4E(I +). (27) 


In order to construct evolutionary tracks in the Hertzsprung-Russell 
diagram it is necessary to evaluate the effective surface temperature Ty, in the 
various cases. Since L~r°T4, this is easily derived. In the early stages of 


evolution 
7 eff 1/2 2 ue 
Tie (2p) (; +X ) 7 (28) 





* The suffix s for surface values will be omitted in this section since M and r always refer to the 
mass and radius of the whole star. 
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After accretion has started to predominate we have 
Ta <A (1+v)"4 
Tog A¥ (2B ong 


The construction of the tracks is now a simple matter. The diagram 
(Fig. 4) represents absolute magnitude against the logarithm of the effective 


temperature. Starting from a definite luminosity L and temperature . at Pijo 
the first slow development forms the straight line Pj). Py. When rapid accretion 
begins, the track turns sharply right. The point of the straight line Pj P4js 
at which this development begins will vary according to the stage of evolution 
along P,). P4), that had been reached when rapid accretion began, 7. e. according 
to the value of yw, Evolutionary tracks PQRS are therefore drawn for 
5 representative values of y, (2, 3, 1, 2, 4). The last forms a limiting case that 
can be almost though not quite reached. For, as was stated before, when X, 
drops below about 0-o1 a considerable part of the energy production will take 
place in the outer zone of the star, and the model here used will no longer be 
applicable. 
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Fic. 4.—Evolutionary tracks in the Hertzsprung—Russell diagram. 
Accreted mass oe 4 : Q 
——————————- = 2, at (2). 
Original mass 1's at R. 


Beyond S part of the outer region is convective. 


As v increases from zero, the representative point of the star will move 
rapidly to the right in the diagram, reaching the neighbourhood of the 
corresponding point Q when v has reached a value of about 0:2. As v increases 
further, the point remains in the same neighbourhood, reaching the point R 
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only when v has reached a value of 1-5. A further increase of v moves the 
representative point upwards and to the left until it reaches the point marked S, 
when a convective zone appears in the outer envelope not far from the 
discontinuity. Although the computations used in this paper no longer apply 
beyond this point, it may easily be seen that in fact the same type of development 
continues. The evolutionary tracks all tend asymptotically to the dotted line 
in the diagram which represents the locus of representative points of 
homogeneous hydrogen stars of differing masses (the main sequence). 

8. The interpretation of Fig. 4 is now quite easy. Assuming a constant rate 
of accretion the representative point would move very rapidly from P to Q, 
and then very slowly in the region Q-R, and would finally move slowly towards S. 
There would hence be a very small density of stars in the P-Q region and a 
high density in the Q-R region. 

If accretion ceases at some stage, mixing becomes again the predominant 
agent of evolution. Since however accretion was supposed to be a much more 
rapid process than mixing, the new developments will proceed much more 
slowly. This argument is reinforced by the consideration that a considerable 
time must elapse before the angular velocity has diffused throughout the star, 
and only then will serious mixing start. The first effect of mixing will be a 
blurring of the discontinuity of composition. Although the effects of such a 
continuous change of composition have not yet been investigated it seems 
unlikely that any substantial change in the structure of the star will occur until 
the composition again approaches uniformity. 

Accordingly the position of the representative point of the star in the 
diagram will be almost stationary after rapid accretion ceases, for a period very 
long compared with the time scale of accretion. 

The developments sketched out presuppose that, while accretion takes 
place, the accretion rate should be high. ‘This in turn implies a low relative 
velocity of star and interstellar cloud, the density of which must not be too small. 
Hence the star is unlikely to move out of the cloud very soon, so that a total 
accretion of less than 20 per cent of the mass is rather unlikely. Since total 
accretions of more than (say) 3 times the original mass are probably not too 
common, we may expect the representative points of stars of the type considered 
to form mainly a broad band along the points marked QR. 

Allowance must now be made for the fact that all stars that underwent such 
a course of evolution did not start off with the same masses. But the chief effect 
of this is merely to expand Fig. 4 vertically. For stars of small mass did not 
have time to convert more than a fraction of their hydrogen into helium so that 
the », cannot greatly exceed 4. Hence at the bottom of the figure only the 
developments corresponding to lower curves can have taken place. Very 
massive stars on the other hand are very uncommon, and in their case high 
helium contents are probably rather likely. Accordingly the upper parts of the 
figure are presumably sparsely populated, with development centring on the 
upper curves. 

Fig. 4 refers, of course, only to stars in which electron scattering predominates, 
but a comparison with Paper II of this series suggests that no very large changes 
are likely to arise if photoelectric opacity is included. For massive red giants 
the low density in any case probably implies that electron scattering predominates. 
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Although further study and detailed comparisons are greatly required, it 
seems eminently likely at this stage that the Hoyle—Lyttleton model of red giants 
offers an adequate explanation of all the main features of that part of the 
Hertzsprung—Russell diagram lying to the right of the main sequence. 

Newnham College, Trinity College, 

Cambridge: Cambridge. 

1951 February 19. 
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ON THE FIGURE OF A PLANET WITH HOMOGENEOUS SHELL 
AND CORE 


Harold Feffreys 
(Received 1951 May 28) 


Summary 


The equations of equilibrium of a two-layer rotating body, with 
uniform densities in the shell and core, are solved without use of the Radau 
approximation. It is found that this approximation loses accuracy 
appreciably when the central condensation is stronger than in the Earth. 
It remains true, however, that the surface ellipticity and the dynamical 
ellipticity are much more sensitive to the ratio C/Ma* than to the radius of 
the core. 





Wiechert’s theory of the figure of the Earth supposed a uniform core 
and a uniform shell, and such models are therefore usually known as of 
Wiechert’s type. His main contribution, however, was to include terms of the 
second order in the ellipticities and to obtain quantitative results in some detail. 
A first-order theory had already been given by Thomson and Tait*, who gave 
the equations but no quantitative results. A still earlier one of S. Haughtont 


is quoted by E. J. Routh}; some quantitative results had been found, the 
author reaching the conclusion, which was considered surprising, that the shell 


must be more than 748 miles thick. 
Such discussions have had little attention since Radau’s discovery that for 


the Earth, on a hydrostatic theory, if 
i ; C-A 


, Fa, (1) 


K = 


~ 2 Ma?’ 
then if m is known, knowledge of any one of K, H and the ellipticity e determines 
the other two within the uncertainty of observation. A model of the type 
considered contains three adjustable constants, restricted only by the need to 
give the right mean density ; then effectively the theory of the figure of the Earth 
gives only one more relation and the problem remains formally insoluble until 
another datum is found. Radau’s theory is a close approximation for a wide 
class of radial distributions of density. It is not universally valid, in fact the 
accuracy begins to fail if K<o-8; and even with the correct value K = , there is 
loss of accuracy if the core radius is less than 2!” of the outer radius.§ ‘Thus 
the approximation might apparently cease to be good with the actual radius of 
the core found from seismology; but Bullen showed || that the variations of the 
density in the core and shell restored its validity. 

* W. Thomson and P. G. Tait, Treatise on Natural Philosophy, Section 831, Cambridge, 1883. 

+ S. Haughton, Trans. Roy. Irish Acad., 22, 251, 1851. 

YE. J. Routh, Analytical Statics, Vol. 2, p. 161, 1902. 

§ H. Jeffreys, M.N., Geophys. Suppl., 1, 121-124, 1924. 

|| K. E. Bullen, M.N., Geophys. Suppl., 3, 401, 1936. 
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It remains interesting, in view of the possibilities of discontinuities of density 
in other planets and stars, to investigate the theory of Wiechert models without 
using the Radau approximation, since some bodies may lie outside its range 
of accuracy. For Jupiter and Saturn, the perturbations of the satellites 
together with the Radau approximation make K near % and } respectively.* 
Using the first-order theoryt and taking the forms of the shell and core 
boundaries as 

r=a(1+eY,); r=aa(1+e'Y,); (2) 
¥,=$— cos* 8, (3) 


and the densities as p, p(I +7), we find, from the conditions that the outer surface 
and the core boundary are level surfaces, 


— fe + na5( —e+ #a7e’) + (1+ na®)m=0, (4) 
Be—(1+4n)e’ +3(1+ ya°)m=0 
pa e+nase’ ee I +9e° 

I+ I+ 


and 


(6) 


For the Earth K is nearly $, and we get, after some algebra, 
ae $a7(I + 2x + 3a7 + a3) 





m  2+42+4%0?+29$03+55}04+ 27905" (7) 
‘The second term is small for all values of « from oto1. This essentially confirms 
Radau’s result from a completely different approach, but our interest is now in 
the amount of the departures. The large coefficients in the denominator are 
much reduced for smaller values of K, and e/m varies more with « for fixed K. We 
also have 

e’ /m = 30a°(1 +«)*/D, (8) 


where D is the denominator in (7). The calculation has been carried out for 
K=}, %, } andi. For K=1, the homogeneous case, e=e’ = H =1-25m for all 
values of «. The results for other values of K are given in Table I in units of the 
fourth decimal. 
TABLE I 
K=3 K=§ 
a e/m fl H/m e/m e’/m H|m e/m e’lm H/m 
°'o 7142 7142 8333 ° 8333 10% ° 104 
orl 7112 35 7040 8306 61 8264 9980 146 9960 
o'2 7022 265 6740 $225 456 8064 9926 1044 9852 
o°%3 6885 798 6283 8085 1434 7751 9860 2786 8 9720 
o"4 6725 1608 5750 7989 2574 7473 9808 4759 9615 
o'5 6578 2437 5198 7886 3880 7219 9776 6422 9553 
06 6444 3475 4774 7810 = 5019-7025 9761 7637 9522 
o"7 6337 4317 4396 7763 5965 6899 9755 8480 9510 
o'8 (6254 5086 4178) 7726 6700 6814 9754 9060 9507 
°'9 (6199 5663 3914) (7705 7263 6763) 9754 9466 9509 
1'o (6154 6154 a (7692 7692 aeac dd (9756 9756 wee 


The shell density is positive when «*<K; values not satisfying this condition 
are indicated by brackets. When «1, with K<1, the shell reduces to a negative 
surface density. 

* H. Jeffreys, M.N., 84, 535, 1924. 
+ H. Jeffreys, The Earth (‘‘ Figure of the Earth”’ chapter), Cambridge, 1924. 


ane 
28% 
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Direct use of the Radau approximation gives 


K=} 5 8 
Na 2°0625 1°25 0°5625 
Wa) O°9171 0°9792 0°9997 
e/m 0°6154 07692 0'9756 
H/m 0°3846 0°6731 0'9512 


As the approximation, in its simplest form, takes ~(y)=1, we should have 
expected that its results would be rough for K=% and seriously wrong 
for K=}. Even for K=§, 7 is monotonic only if «>0-707. In the interval 
0°707 <a <0-913(=K""), where the approximation would be expected to be 
good, it is in fact very close, but outside this interval it is not very good. 

The main interest of this discussion rests in the possibility of determining 
two of the ratios C/Ma?*, e/m and H/m when the other is given; and we find 
that for C/Ma*> a rough determination is possible, but not for smaller 
values. 

An attempt has been made to see whether the present method by itself would 
lead to any useful determination of « when two of the ratios are given. From 
values found in my study* I find for the Earth 

H/m=0-9496 + 0°0002; e/m=0-g9760 + 0°0020. 
The latter depends on the data other than those derived from the constant of 
precession. Taking K=§, a=0°5 as trial values, I find a least squares 
solution 
K =0°8359 + 0°0036; «=0°72 40°13. 

K is nearly right, but « deviates from the true value by more than the standard 
error. This is of course attributable to neglect of second-order terms and of 
differences of density within the shell and core. 


St. Fohn’s College, 
Cambridge : 
1951 May 26. 


*H. Jeffreys, M.N., Geophys. Suppl., §, 219-247, 1948. 





THE SUN’S ROTATION DERIVED FROM SUNSPOTS 1934-1944 
AND ADDITIONAL RESULTS 


H. W. Newton and M. L. Nunn 
(Communicated by the Astronomer Royal) 
(Received 1951 May 26) 


Summary 


Results are given for the Sun’s daily sidereal motion (£) derived from 136 
single recurrent sunspots for the cycle 1934-1944. The expression relating 
€ to solar latitudé (¢) so derived, £=14°-38—2°-96 sin’ ¢, is in satisfactory 
accordance with values from similar data for the five previous 11-year cycles. 
A good internal accordance of measures is found by comparing £ values 
derived from a full rotation period interval (27 days) with the mean (£,), given 
by the two separate disk-passages of the same spot. 

From the combined data of six sunspot cycles there appears no 
measurable variation in the Sun’s rotation period (1) with alternate 11-year 
cycles, (2) within the (mean) cycle itself. Neither do the two hemispheres, 
north and south, behave differently. 

Although not included in the discussion, the derived expression for the 
daily sidereal motion from non-recurrent spots, 1934-1944, is stated as a 
point of interest in the Addendum. 





Introduction.—The problem of the Sun’s rotation still deserves more 
observational data. Considerable differences in spectroscopic values cannot 
wholly be ascribed to differences of level in the chromosphere for the particular 
line or lines used. Moreover, St. John pointed out in 1932* that collective 
spectroscopic results ‘‘show a progressive change in the equatorial rotation 
determined from the reversing layer during the past thirty years greater than 
can be attributed to errors of observation’’. The overall range in the (smoothed) 
observational values from 1906 to about 1925 is equivalent to no less than two 
days in the respective rotation periods for the Sun’s equator. 

Values derived from sunspots are more accordant inter se. The range in the 
equatorial value (a) for the daily sidereal rotation (é) at the equator from the 
expression £=a—bsin*¢ given by some half-dozen investigations ¢ is 0°22, 
equivalent to about 0-4 day in the rotation period. 

Very consistent values over five sunspot cycles (1878-1933) have been given 
by the Greenwich measures of individual spots of long duration.{ The range 
in the value of a does not exceed 0°-04 per day. The dispersions of € values in 
the 5° zones of latitude, into which the data are divided, is also much less than 
in determinations based on spot groups (as in the investigation by Mr and 
Mrs Maunder §), especially those that are shorter-lived. 

In the present paper, the results from recurrent spots for a sixth sunspot 
cycle (1934-1944) are presented. Over the whole period since 1878 there are 
now sufficient homogeneous data from at least two of the 5° latitude zones to 

*C. E. St. John, Trans, I.A.U., 4, 43, 1933; see also P.A.S.P., 47, 295, 1935. 
t The smallest value of a is 14°-35 (Spoerer) and the largest 14°57 (O. A. Akesson). 
t Royal Observatory, Greenwich, M.N., 85, 548, 1925. H.W. Newton, M.N., 95, 62, 1934. 


§E. W. and A. S. D. Maunder, M.N., 65, 813, 1905. O. A. Akesson, Medd. Lunds Astro- 
nomiska Observatorium, Serie II, Nr. 11, 1913. 
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show whether £ varies within measurable limits (1) alternately from one cycle 
to the next, as do the magnetic polarities of sunspots, and (2) with the phase of 
the (mean) 11-year sunspot cycle. Results of this enquiry are stated herein. 

Data.—The basic data comprise individual recurrent sunspots tabulated in 
Ledger I of the Greenwich Photoheliographic Results for the years 1934 to 1944. 
The volumes for 1937 and 1939 are in course of publication, those for 1940 to 
1944 being still in MS. 

Although the advantages in using stable recurrent spots appear far to outweigh 
the disadvantages, it is true nevertheless that the non-uniform characteristic 
motions*, especially of the leader spots of groups, are imposed on the true values 
of £. The present data include a large proportion of leader spots first seen as 
such or later as the sole survivor of the group, when they would be classified as 
“‘unipolar’’ spots. 

Counting each appearance separately, the spots used in this paper were as 
follows :— 


Leader spots (89): Followers (8): Other components (1): 
Single spots, i.e. ‘‘ unipolar’’ (174). 

Any systematic effect from non-uniform proper motion on values of € is 
much reduced by (1) the data including spots in all stages of development, and 
(2) the use of the long time-interval of a solar rotation for deriving the values of €. 

Regarding the stage of spot development, the distribution in three broad 
divisions for the present data is as follows :— 

Number Mean areat 
Phase I Origin or early development to maximum 63 394 
II Middle life 139 298 
III Late decline (and extinction) 7O 71 


The time spent in each of the three phases is very approximately in the ratio I: 2: 1. 
Phase II is sometimes markedly prolonged with a spot of area 200 or 300 millionths. 

Derivation of £.—The method of reduction as used in earlier Greenwich papers 
is simply as follows :— 

The mean longitude of the spot during its first complete (or part) disk-passage 
is compared with the corresponding mean longitude at its second appearance. 
Days on which any spot was >80° from the central meridian are excluded. 

If the change in longitude be x° in the interval y mean solar days, then the 
diurnal angular motion is x°/y relative to a meridian on the Sun having a sidereal 
rotation of 25-38 days } or 14°18, per day, this being the datum adopted in the 
Nautical Almanac for computing the daily values of the longitude of the central 
meridian. ‘Thus the daily angular sidereal motion, £, per mean solar day is 
14°°18+x°/y. In rounding off mean values of € to two decimal places of a 
degree (as in Table I), the correction of +0°-004, being the excess over 14°18, 
has been taken into account. 

Observational checks.—The comparison of the mean position of the same 
sunspot 27 or 28 days apart implies a (close) symmetrical distribution of the 
measures of longitude with respect to the central meridian. In the case of 

* Royal Observatory, Greenwich, M.N., 85, 553, 1925. 
+ Areas, corrected for foreshortening, are expressed in millionths of the Sun’s hemisphere. 
} It is perhaps worth recording that this was the period adopted by Carrington because of “‘ its 


admitting conveniently of much subdivision without remainders ’’’. He found later that this value 
applied to latitude 14°. 
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incomplete disk-passages of spots, the daily longitudes for comparison at the 
two epochs were selected so as to make the mean distances from the central 
meridian (and on the same side of it) correspond as nearly as possible. In this 
way, any small systematic displacements that may be present in the measures 
and which are a function of the distance from the centre of the disk will be 
eliminated. A small systematic error of this nature—zero at the centre of the 
Sun’s disk and a maximum at the limb—could arise from differential fore- 
shortening (not allowed for by the measurer) of the area subtended by the larger 
spots approaching the limb. This will tend to make the measured distance, r, 
of the spot from the disk centre smaller than its true value. 

Other displacements have been suggested * as arising from (a) refraction 
above the photosphere producing an augmentation of the Sun’s radius measured 
as R (irradiation will also produce an augmentation), and (4) a possible ‘ depth- 
parallax”’ arising from a difference of level between spots and the effective limb. 
The combined effect of these possible displacements would tend to decrease the 
ratio r/R which enters directly into the computation of longitude and latitude 
of the spot, and thus decreases the derived daily angular motion. Although the 
use of recurrent spots avoids these systematic errors inherent in single disk- 
passage determinations, there is occasional uncertainty in the continuity of the 
longer-lived spots from one rotation to another. A general check is therefore 
desirable by comparing € values given by the same spot from the 27-day interval 
and the mean value (é,) given by the separate disk-passages. Cases of complete 
or nearly complete disk-passages only have been included for this purpose. 
For each of the two disk-passages the mean longitude from the first five days is 
compared with that from the last five days. Days on which the spot position 
exceeded 80° from the central meridian are excluded, so that the reduced effect 
of the small possible displacements considered above would be unlikely to 
invalidate the check, as such, on the continuity of any particular sunspot. 

The distribution of values, £—£,, from 86 pairs of disk-passages is shown in 
Fig. 1, which also incorporates the development phases of the spots concerned. 
With respect to the centre grouping of values, say from —0°-10 to +0°*I0 
(comprising 63 of the 86 values), it may be said that the distribution, when 
slightly smoothed, is essentially that of a normal error distribution for a probable 
error of + 0°-03 (single value) and displaced +0°-02 from the origin. 

This small positive asymmetry of 0°-02, which is five times the p.e. of the 
mean of the above 63 observations, may be real and represent the combined 
effects of the small systematic errors entering into a single disk-passage 
determination of the solar rotation. 

It may be added that differential foreshortening alone would amount, on 
purely geometrical grounds, to about 0°-or in the above derivation of £, values 
from circular spots of the average order of 3° in diameter or 350 millionths of 
the Sun’s hemisphere in area.t A systematic error of this order could well 
arise from the optical distortion of the lens system of a photoheliograph, as 
instanced by the small corrections used since 1926 on this account for the 
Greenwich photographs. 


* Papers of the International Union of Solar Research, Computing Bureau: R. S. Capon, 
M.N., 73, 361 and 732, 1913, in which various references are included. 

t We are indebted to Mr H. F. Finch for the evaluation of this effect after an informal discussion 
on possible systematic errors. 
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The large negative values in Fig. 1 are found to be all* associated with 
Phase I of leader spots, when their forward motion in longitude would entirely 
obscure the other very small displacements. 

A few other high values of £—£, indicated doubtful continuity of the spots 
on the Sun’s invisible hemisphere, and this uncertainty was confirmed by a 
re-examination of their daily positions and areas. Such cases are not included 
in Fig. nor in the data from which ¢ was derived. 


+ —°20 —-10 0 +10 +20 


(> 20) 
Fic. 1.—Frequency distribution of §—£, values and phase development of the associated recurrent 

sunspots at both disk-passages. 
Ordinates : Frequency. 
Abscissae: Values of §—€,. 
£] Phase I. 
—} Phase II. 
fA Phase III. 


The overall mean, £ —£,, is less than 0°-or (—0°-007). This small statistical 
difference not only substantiates the continuity of the spots provided by Ledger I, 
but suggests that non-recurrent spots from Ledger II can serve at least as a 
useful independent check for Table I. 

The 282 non-recurrent spots selected include some from Ledger I which do 
not meet the requirements of a 27-day interval comparison. The derived values 
of the solar diurnal motion from these single disk-passages (designated £&,) are 
tabulated in the Addendum. 

Results,—In Fig. 2 the separate values of € from 136 recurrent sunspots are 
plotted against their respective solar latitudes (¢). ‘The same values averaged 
for 5° zones of latitude (0°o to 4°-9; 5°:0 to g°-9, etc.) are given in Table I with 
northern and southern hemispheres shown separately in satisfactory accordance. 

* Each recurrent spot is of course assigned two indices of phase development. Phase I scarcely 


ever extends to the second passage. Hence in Fig. 1, a 50 per cent occurrence of Phase I indicates 
a complete association of these spots with Phase I at their first disk-passage. 
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SS ae eee s re es = i i l 
0° 5° 20° a 30° 35° 
Fic. 2.—Observed values of & from recurrent sunspots, 1934-1944. Means for each 5° zone of 
latitude are represented by open circles. The curve is computed from the expression 
&= 14°°38—2°-96 sin? ¢. 
Ordinates (left-hand side): Daily sidereal motion. 
Ordinates (right-hand side): Sidereal rotation period in days. 
Abscissae: Solar latitude. 
Tas.e I 
Mean daily sidereal motions of sunspots, 1934-1944 
Mean 
N. spots 5S. spots Mean dispersion O-C 
3°°4 4°°0 3°°6 
14°°40 14°°43 14°-41 + 0°-02 + 0°08 +0°04 
6 16 
7°°8 7°°8 
14°34 =. 14°*33.0°-01 + 0°09 
14 29 
11°'9 12°°4 
14°°23 14°:22+0°-o1 + 0°08 
18 36 
17°°2 197°°2 
14°12 14°-12+0°:o1 +0°-06 
15 22 





= moe 2 Ho 


wn 


~ 
= 


22°°7 22°°3 
13°90 13°*96+ 0°02 + 0°08 
9 20 
27°°3 27°°6 
13°°77 13°77 + 0°03 +0°*10 
8 II 
32°°4 
13°48 +0°14 
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In column 5 the mean dispersion (or scatter of observed points in Fig. 1) is given, 
from which the respective probable errors in the preceding column have been 
computed. 

Taking the observed mean values of £ in each 5° belt of latitude and weighting 
the values according to the respective number of observations (the mean 
dispersion in each 5° belt being comparable), a solution by least squares was 
obtained, giving 

£ =14°-38, + 0°-009 — (2°96 + 0°-09) sin® d. 


Residuals (observed minus computed) are given in the last column of the table. 
This evaluation of € is in good accordance with those for the five earlier sunspot 
cycles as shown in Table II. 
In Table II this result is compared with those yielded by the preceding five 
cycles.* 
Taste II 
Solutions for & (recurrent sunspot data) over six cycles 


Number of Diurnal motion (&) 
Cycle Years sunspots 
I 1878-1888 89 14°°34+ 0°01 —(2°'35 + 0°08) sin? 
II 1889-1899 122 14°°36+ 0°01 —(2°*58+ 0°18) sin? d 
Ill 1900-1913 99 14°°38+ 0°01 —(2°-69+ 0°13) sin? d 
IV 1914-1923 139 14°°38+ 0°01 —(2°-50+ 0°18) sin? d 
V 1924-1933 130 14°°37+0°-01 —(3°°01 +011) sin?’ d 
Vv 1934-1944 136 14°°38+ 0°01 —(2°-96+ 0°09) sin? 


Using the present solution (6th cycle), values of € at 5° steps of latitude are 
given in Table III, together with the equivalent rotation periods. These values 
will be found to compare best with those derived by Carrington from his visual 
measures of sunspots 1853 to 1861. 


TasLe III 
Computed values of £ and equivalent 
rotation periods from = 14°:38,—2°'96 sin? ¢ 


Solar Daily sidereal Sidereal Mean synodic 
latitude motion, £ period period 
d 


° d 
14°38 25°03 26°87 
5 14°36 25°07 26°92 


° 


10° 14°29 25°19 27°06 
14°18 25°38 27°28 
20° 14°04 25°65 27°59 
25° 13°85 25°99 27°98 
30° 13°64 26°39 28°45 
35° 13°41 26°85 28°98 


Test for periodic variation of £.—(1) Alternation with cycle; (2) during mean. 

cycle. 
(1) Alternation with cycle.-—The mean values of £ for the odd and even cycles 
have been taken from Table II and weighted according to the respective relative 
* These solutions differ slightly from those given in M.N., 95, 62, 1934 (Table II) for which 


a less rigorous solution was used. In that table, lines 1 and 2, read 14°35 and 14°38 respectively 
for 14°°36 and 14°39. 
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weights of a and 4 in each of the six solutions. The derived expressions, in close 
accordance, are as follows :— 

Cycles I, III, V. £ = 14°36, + 0°-005 — (2°60 + 0°06) sin? ¢, 

Cycles 11, IV, VI. €=14°-37;+0°-007 —(2°-82 + 0°07) sin® d. 

However, the derivation of the equatorial value from an expression which 
best fits the whole range of latitude values (with considerable weight because 
of numbers in the middle latitudes) might obscure some small significant 
differences near the equator itself. We therefore give the observed values of & 
for the equatorial zone, 0°—5°, and also these values reduced to the equator by 
the present solution for €. 

Observed & d £ at o° 
Odd cycles 14°°369+ 0° ols 3°°3 14°36, 
Even cycles 14°°40,+ 0°01, 2°"9 14°°41, 
The difference of 0°-05, being less than three times its probable error, is therefore 
of very doubtful significance. 

The two most prolific 5° zones of latitude (10°-15°; 15°-20°) provide the 
following observed mean values over the same extended period, 1878-1944. 
Values of € in italics indicate the slightly adjusted values to a common mean 
latitude of 12°-5 for the lower zone and to 17°-2 for the higher. m is the number 
of spots. 

TABLE IV 
Cycles Zone 10°—15° Zone 15°-20° 


é $4 : $ on 

Odd 14° 236 12°°7 86 14°°13, s7""@. 93 
14°"249 14°*I24 

Even 14°'239 sa°3 42g 14°*I2, 17 
14°22, 14°*I39 


o 


%:- @B 


From these values we see that there is no observational evidence ‘of any 
cyclical variation of £ in these latitude ranges as great as I part in 1000 in a 22-year 
period. 

(2) Variation of & within the 11-year cycle.—To test whether there is any 
periodic variation of € during the (mean) sunspot cycle, the data for the above 
two zones were divided into four epochs according to time-interval from solar 
maximum. ‘The epochs of the six maxima were adopted from the Ziirich data 
of sunspot “‘ numbers ”’, and the limits of each of the four epochs within the mean 
11-year cycle were adjusted so as to include as nearly as possible the same number 
of observations in the same 5° zone. ‘The results are set forth in Table V. For 
more exact comparison, values of € in italics include the small correction necessary 
to bring the values appropriate to the mean latitude of 12°-5 in the case of the 
first zone and of 17°-2 for the other zone. Although the greatest difference 
between one epoch and the next is 0°-08 (i.e. about four times the probable 
error), the overall run of values suggests no real variation with phase of the cycle. 

This result from sunspots (with the reservation that the solar minimum itself 
is barely represented) contrasts sharply with those typical of the upper chromo- 
sphere for which Evershed*, for instance, obtained from measures of the H and 
K lines in prominence spectra, a difference of 2° per day as between solar 
minimum and times of greatest activity. 


* J. Evershed, M.N., 105, 205, 1945. 
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TABLE V 


€ values during the (mean) sunspot cycle 
Epoch from sunspot maximum (=o0*-o) 


I 2 3 
—19-7 —oy'1 +193 
12°°6 12°°5 12°°5 
14°'22 14°°26 14°°22 
14°°22 14°°26 14°22 
52 51 5° 


—29+3 —1Y'o 4 oy's 
197°°9 17°°8 16°°8 
14°°17 14°11 14°°1I 
14°°18 14°°13 14°°10 

n 34 37 37 


Mean é, (10°-15°) 14°'233 at latitude 12°-5 (201 spots). 
1878-1944 | (15°-20°) 14°-126 at latitude 17°-2 (146 spots). 


Northern and southern hemispheres compared (1878-1944).—Finally from the 
combined data, a solution for € has been made for the separate hemispheres * 
to see whether there is any measurable difference in behaviour. It will be seen 
that there is no difference within the range of present accuracy. 


North: €=14°-37, + 0°-010 —(2°-68 +.0°-12) sin? (357 spots), 
South: £=14°-38, + 0°-005 —(2°-91 +0°-07) sin? (358 spots). 

Combining data for six cycles by 5° latitude zones and re-solving, we obtain 

£=14°-37, + 0°-006 — (2°77 + 0°-08) sin? ¢. 

It is of interest to compare the above solution with that obtained from the 
six individual solutions in Table II combined according to the weights of a and b, 
as follows :— 

& = 14°36, + 0°-004 — (2°69 + 0°-04) sin? d. 


Acknowledgment is gratefully made to the Astronomer Royal for the 
opportunity given to prepare this paper; for advice concerning the weighted 
solution adopted, and also for his suggestion that the data might be arranged to 
show any existing periodic variation of £. 

With the kind cooperation of the Superintendent, H.M. Nautical Almanac 
Office, the various data on the original working cards were transferred to punched 
cards and so passed through the Holleritn machine to give the collected values 
in Table I and in the Addendum, together with other information. 


Addendum 


Using data from the non-recurrent sunspots for the period 1934-1944, 
values (designated €,) of the daily sidereal motion have been derived from 5-day 
means of longitude as for ¢, values (i. e. separate disk-passages of recurrent spots). 

Of the 282 spots available, mainly from Ledger II of the Greenwich Photo- 
heliographic Results but including a few components of groups in Ledger I that 
were not suitable for the basic discussion, 89 spots were seen as leaders of groups, 


* As for the first four cycles in M.N., 85, 552, 1925. 
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41 as followers, three as other components and 149 as single or “ unipolar”’ spots. 


The distribution of their development-phases is as follows, the mean areas being 
added in brackets :— 


Phase I—20 (309); Phase II—170 (236); Phase III—g2 (122). 
The above phase assessment is necessarily more approximate than for the 
long-lived spots. Sometimes, the whole life-history of the spot was completed 
during its disk-passage, and in such cases the index II was assigned. 


Mean values of £, in 5° belts of latitude were derived as for € values in Table I, 
and a solution made as before, which gave :— 


£ =14°-41, + 0°-004 — (3°36 + 0°-04) sin’ d. 
The residuals (O—C) for the respective 5° latitude belts are satisfactorily small 


and read as follows from the equatorial belt to that for 30°-35°, the respective 
numbers of spots being added in brackets :— 
+0°02 (9); +0°o1 (58); —o°or (87); 0°00 (54); —o°or (41); 
+0°-02 (20); and 0°00 (13). 

It seems worth recording the above subsidiary result without stressing the 
possible significance of the larger values of the a and 5 coefficients as compared 
with those derived from long-lived spots. Although a small physical difference 
(e. g. level difference) between these two classifications of spots may be indicated, 
it is nevertheless true that the values derived for the solar rotation depend to 
some extent on the selection of spots used and their phase and rate of development. 
A more rapid development of the shorter-lived spots may contribute larger 


values of £, (contrary to the effect of small systematic errors previously mentioned 
on p. 415), since the forward movements in longitude of leaders during growth 
will be effectively averaged over a 9 or 10 day interval instead of 27 days, as 
for the recurrent spots. 


Royal Greenwich Observatory, 
Herstmonceux Castle, 
Sussex : 
1951 May 19. 
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The following results are in continuation of those given in M.N., 110, 501, 
1950, and are derived from the measurement at Greenwich and Herstmonceux 
of photographs taken at the Royal Observatories of Greenwich and the Cape and 
at the Kodaikanal Observatory, India. 

Table I gives the mean daily areas of umbrae, whole spots and faculae for 
each synodic rotation of the Sun included in the year 1946; the means for each 
year as a whole are included in Table II, which summarizes the yearly values 
since the last sunspot minimum of 1944. 


TABLE I 
Mean Daily Areas 


Projected* Corrected for Foreshorteningt 
No. Rotation Days. ———____~—- a polation VE carreras 
of Commenced Photo- Whole i Whole 
meontiats U.T. graphed Umbrae Spots Faculae Umbrae Spots Faculae 
1946 
1235 Jan. 4°63 27 295 1648 1093 242 1376 1353 
1236 §©Jan. 31°97 27 709 3886 1619 550 3108 1917 
1237 ‘eb. 28°31 28 306 1882 1846 255 1600 2210 
1238 Mar. 27°62 27 296 1695 1826 240 1401 2201 
1239 =. Apr. 23°89 27 299 1806 1699 210 1282 2074 
1240 May 21‘12 27 149 854 1653 III 649 2049 
1241 June 17°32 28 247 1526 1844 184 1143 2285 
1242 July 14°52 27 727 5170 =: 1905 537 3915 2427 
1243 Aug. 10°74 27 439 2534 1600 311 1836 1966 
1244 Sept. 6°98 27 377 2351 1799 276 1732 2185 
1245 Oct. 4°26 28 328 1924 2065 236 1387 2552 
1246 Oct. 31°55 27 418 2629 2046 295 1896 2489 
1247. Nov. 27°86 27 489 3167 2338 359 2392 2815 
* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 


TaBLe II 
Mean Areas and Heliographic Latitudes 
Mean Daily Areas 


No. of Days Projected* Corrected for Foreshorteningt 


Umbrae Whole Faculae Umbrae Whole Faculae 


Photo- Without 
graphed Spots Spots Spots 
366 157 30 160 281 23 126 338 
365 14 102 560 773 78 429 939 
365 ° 389 2381 1793 291 1817 2188 


* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 
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Table III gives for each rotation in the year 1946 the mean daily area of the 
whole spots (corrected for foreshortening) and the mean heliographic latitude of 
the spotted areas for both northern and southern hemispheres. 

The mean heliographic latitude of theentire spotted area andthe mean distance 
from the equator of all spots are also tabulated. ‘The mean values for the year 
1946 are included in Table IV. 

Tables II and IV are in continuation of similar tables in Monthly Notices: 
for the years 1874 to 1888, 49, 381, 1889: 1889 to 1902, 63, 465, 1903: IgOr to 
1914, 76, 402, 1916: 1913 to 1924, 85, 1007, 1925: 1923 to 1933, 94, 870, 1934; 
and 1933 to 1945, 110, 50I, 1950. 

The rotations in Tables I and III are numbered in continuation of Carrington’s 
series, 

The chief features of the record are as follows :— 


(1) The increasing steepness in the rise in sunspot activity after a solar 
minimum in 1944:2. As compared with the year 1945, the mean daily sunspot 
area has risen more than four times and the faculae 2} times (Table II). 

The large mean daily sunspot area value of 3108 millionths (corrected area) 
of rotation 1236 (Table I) is largely due to the giant sunspot which, although 
supplanted in 1947, was up to this date the largest ever recorded at Greenwich. 
Similarly, a second giant spot in July contributed largely to the high value of 
3915 in rotation 1242. 

The run of the mean daily areas of faculae from rotation to rotation is smoother, 
because of the much longer life of faculae arising with spots. 

TABLE III 
Spots North Spots South 
of the Equator of the Equator Mean 
——~ Mean Distance 


Mean Me: Mean Latitude from 
Helio- “*©8" Helio- of Entire Equator 


No. Rotation Mean 
of Commenced Daily 


aily 
Rotation ti Area 


raphi 
ee a Area 


Latitude 


graphic Spotted of all 
Latitude Area Spots 


1946 : p 
1235 Jan. 4°63 905 22°15 471 16:29 §=6©++ 8-99 20°14 
1236 Jan. 31°97 2514 25°44 594 19°98 = + 16°77 24°40 
1237 Feb. 28-31 943 25°72 657 22°73 + 5°81 24°49 
1238 Mar. 27°62 510 26°60 891 30°48 — 97! 29°07 
1239 Apr. 23°89 487 24°58 796 21°54 — 4°03 22°69 
1240 May 21°12 510 19°83 139 20°30 +-11°26 19°93 
1241 June 17°32 522 24°10 620 22°81 - 1°37 23°40 
1242 July 14°52 2721 20°24 1194 21°89 + 7°39 20°74 
1243 Aug. 10°74 1434 18°04 401 16°70 + 10°45 17°75 
1244 Sept. 6°98 1295 18°17 436 15 ‘04 + 9°81 17°38 
1245 Oct. 4°26 519 15°69 868 14°83 — 3°41 15°15 
1246 Oct. 31°55 1192 15°27 704 13°31 + 4°66 14°54 
1247 Nov. 27°86 1052 14°93 1341 10°31 + 0°78 12°34 

(2) On no day in 1946 were sunspots and faculae absent from the Sun’s disk. 
(3) The mean daily area of the northern hemisphere was greatly in excess of 
that of the southern. But of the 437 units of the northern excess, 309 units 
were alone due to the two giant spots mentioned above. In actual numbers of 


groups, the southern was rather the more prolific, as shown next. 
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TABLE IV 
Spots North Spots South 
of the Equator of the Equator Mean 
a i Mean Distance 
Mean Mean Latitude from 
Helio- Mea 2 Helio- of Entire Equator 
graphic Daily graphic Spotted of all 
Latitude AT Latitude Area Spots 


Mean 
Daily 
Area 


Days 
Photographed 


° 


1944 366 42 19°00 83 22°81 —8-70 21°53 


Old Cycle 7 4°17 6 7°61 —I'Il4 5°72 
1944 366 New Cycle 35 22°18 77 24°02 —9°63 23°45 
1945 365 121 20°13 309 20:26 —8-92 20°22 
1946 365 1127 20°74 690 18-79 +5°73 20°00 


(4) The number and distribution, northern and southern hemisphere, of spot 
groups of (a) two days or longer; (5) one day only, are as follows :— 


(a) (5) 
Northern spots 179 32 
Southern spots 190 43 


Total 369 75 


The range in latitude for the longer-lived spots was from 6° to 42°. The latitudes 
of the 1-day spots ranged from 1° to 41°. 

(5) The mean weighted latitude of all spots was 20°-o—practically the same 
as for 1945. ‘The latitudes of the two giant spots were +26°-I and +22°-2 
respectively. 

A short descriptive survey of the sunspots of 1946 has been given in Monthly 
Notices, 107, 103, 1947. 

An appended table gives the mean daily areas of sunspots and of faculae 
(corrected for foreshortening and expressed in millionths of the Sun’s hemisphere) 
for each calendar month of 1946. 


Monthly Mean Daily Areas of Sunspots and Faculae 
1946 

Month Spots Faculae Month Spots Faculae 
Jan. 1237 1285 July 2912 2431 

Feb. 3091 1966 Aug. 2291 2010 

Mar. 1551 2283 Sept. 1860 2181 

Apr. 1433 2098 Oct. 1345 2564 

May 1085 2036 Nov. 1852 2458 

June 716 2099 Dec. 2514 2819 


Royal Greenwich Observatory, 
Herstmonceux Castle, Sussex : 
1951 March. 
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